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ABSTRACT 

We present a detailed assessment of the various dynamical pathways leading to the coalescence of 
compact objects in Globular Clusters (GCs) and Short Gamma-Ray Burst (SGRB) production. We 
consider primordial binaries, dynamically formed binaries (through tidal two-body and three-body 
exchange interactions) and direct impacts of compact objects (WD/NS/BH). Here we show that if 
the primordial binary fraction is small, close encounters dominate the production rate of coalescing 
compact systems. We find that the two dominant channels are the interaction of field NSs with 
dynamically formed binaries, and two-body encounters. Under such conditions, we estimate the 
rcdshift distribution and host galaxy demographics of SGRB progenitors, and find that GCs can 
provide a significant contribution to the overall observed rate. 

Regarding the newly identified channel of close stellar encounters involving WD/NS/BH, we have 
carried out precise modeling of the hydrodynamical evolution, giving us a detailed description of the 
resulting merged system. Our calculations show that there is in principle no problem in accounting for 
the global energy budget of a typical SGRB. The particulars of each encounter, however, are variable 
in several aspects, and can lead to interesting diversity. First and most importantly, the characteristics 
of the encounter are highly dependent on the impact parameter. This is in contrast to the merger 
scenario, where the masses of the compact objects dictate a typical length and luminosity scale for 
SGRB activity. Second, the nature of the compact star itself can produce very different outcomes. 
Finally, the presence of tidal tails in which material will fall back onto the central object at a later 
time is a robust feature of the present set of calculations. The mass involved in these structures is 
considerably larger than for binary mergers. It is thus possible to account generically in this scenario 
for a prompt episode of energy release, as well as for activity many dynamical time scales later. 
Subject headings: gamma rays: bursts — stars: neutron — hydrodynamics — accretion disks - 
gravitational waves — globular clusters 



1. INTRODUCTION 

Until fairly recently, SGRBs were known predomi- 
nantly as bursts of 7-rays, and largely devoid of ob- 
servable traces at lower energies. The launch and 
successful operation of the Swift satellite has now en- 
abled the detection and localization of X-ray after- 
glows from several events, enabling in turn the study 
of their properties at optical and radio wavelengths 
and the identifica tion of the host g a laxies at cosmo - 
logical distances (IBerger et alJ 120051 : IFox et all I200E 



Barth elmv et alj|2005t iHiorth et all 12005] IGehrels et alJ 
20051: iBloom et al.ll2006HProchaska et al.l(2006t) . The oc- 



currence of a fraction of events among old stellar popu- 
lations e.g., of an elliptical galaxy for GRB 050724, rules 
out a source uniquely associated with recent star for- 
mation. In addition, no bright supernova is observed 
to accompany SGRBs, contrary to what is seen in most 
nearby long-duration GRBs. It is clear by now that short 
and long events are not drawn from the same parent stel- 
lar population, and that SGRBs are far from st andard 
(jNakarl 120071 : IGehrels. Ramirez-Ruiz fc Foxl[200l . This 
hints at the underlying possibility that the progenitor it- 
self may be quite different from burst to burst, and not 
entirely restricted to the most widely favour ed scenario 
involving the merger of close binaries (e.g. iPaczvriskl 



H986tlEichlerlll989l ) containing neutron stars (NS) and/or 
black holes (BH). 

Since the first evidence from the Uhuru and OSO-7 
satellites revealed a population of highly luminous low- 
mass X-ray binaries (LMXBs) in globular clusters (GCs), 
it has been noted that the formation rate per unit mass 
of these objects is orders of magnitude higher i n GCs 



ot these objects is orders ot magnitude higher m UUs 
than in the Galactic disk (|Katzlll975t lClarkj|l975l) . This 
discovery stimulated a flurry of theoretical work into the 
formation of GCs LM XBs by the process e s of two- and 
three-body encounters (|Fabian et al.l Il975 : Heggic 1975; 
lHillsl ll975). These dynamical formation scenarios are a 
natural explanation for the high occurrence of LMXBs 
in GCs since the stellar densities, and hence encounter 
rates, are much higher in the cores of GCs than other 
regio ns of the Galaxy (jPoolev et alJ 120031 : iHeinke et alJ 
2006). In such environments, in fact, it is unavoidable 
that many stars undergo close encounters and even physi- 
cal collisi ons, with high probability, within their lif etimes 
(see e.g., iRosswog et ail 120091: iRaskin et ai]|2009D . It is 
the interplay between compact stars in such dense envi- 
ronments and their ability to trigger SGRBs that forms 
the main topic of this paper. 

In general, forming a compact binary system requires 
some mechanism to dispose of enough energy to effec- 
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tively bind them. Two such processes h ave been consid- 
ered for conditions in globular clusters (|Hut fc Verbunti 
1983). The first is the presence of an additional star 
which can carry away s ome energy in k i netic form. This 
has been considered by iGrindlav et al.1 ((2006) as a plau- 
sible SGRB channel by computing the outcomes of ex- 
change interactions between binaries containing one com- 
pact object and a single neutron star in a collapsed GC 
core. The second is the loss of energy to internal stellar 
oscillations, excited by the tidal forces of one star on the 
other. 

Here we suggest an alternative and perhaps less re- 
strictive mechanism for SGRB production related to the 
second mechanism given above, namely, the tidal cap- 
ture and collision of compact objects in dense stellar en- 
vironments. In this new scenario, the compact objects 
are contained within a globular cluster, and interact di- 
rectly through close encounters rather than being driven 
together by pure gravitational wave emission in existing 
close binaries. Event rates for such interactions within 
galaxie s have been found to b e much too low to be of in- 
terest (jJanka fc Ruf fert 1996|) when compared with the 
GRB event rate, but they may be frequent enough to 
have an important effect on their production in GCfQ. 
In this paper we make a careful assessment of the vari- 
ous dynamical evolutionary pathways involving compact 
objects in such environments, focusing particularly on 
the frequency and physical character of the tidal cap- 
ture of two compact objects that are ultimately capable 
of leading to a merger and thus powering a gamma-ray 
burst. We find that they can provide a substantial con- 
tribution to the total rate, and compute the evolution of 
the associated rates with redshift. 

Some pressing questions include: When would two 
passing relativistic stars capture each other in a bound 
close orbit? When two neutron stars collide, does the 
rapidly-spinning merged system have too much mass (for 
most presumed equations of state) to form a single sta- 
ble object? If so, the expected outcome after a few mil- 
liseconds would therefore be a spinning BH, orbited by 
a torus of neutron-density matter. When a NS collides 
with a BH, does enough mass remain in the orbiting de- 
bris to catalyze the extraction of energy from the hole 
at a rate adequate to power a short-lived GRB? How 
do the long tidal tails thrown out through the outer La- 
grange point affect the accretion stream around the pri- 
mary star? Even if the evolution time scale for the bulk of 
the debris torus were no more than a second, is enough 
mass and energy still available to power the late time 
flares? What is the relative frequency and observable 
signatures of all these collision events and how do they 
compare to Swift observations? 

The structure of this paper is as follows: detailed hy- 
drodynamic simulations of encounters of compact objects 
of various types and with varying impact parameters are 
presented in Section [5] together with a detailed descrip- 
tion of the numerical methods and the initial models; 
the resulting gravitational wave signals are shown in Sec- 
tion [3] Section |4] offers an estimate of the encounter rate 

1 Hansen & Murali (1998) have proposed collisions of compact 
objects with main sequence stars as possibly relevant for the pro- 
duction of GRBs in GCs, although it is not clear how the large 
baryon loading problem could be circumvented in this case 



as function of cosmic time, and compares it with the 
merger rate of compact binaries. We discuss our find- 
ings in Section [5j and we summarize and conclude in 
Section [6] 

2. DYNAMICS OF TIDAL CAPTURE AND DISRUPTION 

2.1. Numerical implementation 

The hydrodynamical calculations shown subsequently 
have been performed with the three dimensional Smooth 
Particle Hydrodynamics (SPH) code previously used 
to study merging dou ble neutron star and black hole- 
neutron star binaries (|Monaghanlll992t iLee fc Kluzniakl 
Il999h . The problem does not allow for simplifications due 
to symmetry and SPH, being a Lagrangian scheme lends 
itself particularly well to this kind of sit uation. The tidal 
tails s e en in bina r y mergers [see e .g. | Rasio fc Shapirol 
(|1994D ;lLei (|200lD : iRosswog et all (|2003D ] are a natural 
outcome in the present scenario as well, and following 
their formation and evolution is one of the main objec- 
tives here, which is not possible with grid-based codes. 
Self-gravitating spherical stars of a given mass, M, and 
radius, R, in hydrostatic equilibrium with a polytropic 
pressure-density relation P = K p T , where K and T are 
constants, are constructed with N ~ 10 5 fluid elements 
(SPH "particles") and used for the dynamical evolution. 
We have performed test simulations with varying spatial 
resolution to test for convergence, using from N ~ 10 4 
to N ~ 2 x 10 5 SPH particles. We are here mainly in- 
terested in the dynamics of the disruption process for a 
limited number of dynamical times, and have found that 
N ~ 10 5 particles are sufficient for this purpose. We 
thus report all simulations at this initial uniform res- 
olution level. Our calculations are Newtonian, so we 
cannot model a true black hole. We merely approxi- 
mate one as a point mass Mbh, any matter approach- 
ing within a Schwarzschild radius rg c h = 2GMbh/c 2 be- 
ing accreted (and its mass added to that of the hole). 
Clearly this is a problem that ideally should be treated 
in full General Relativity, but a first understanding can 
be gained with the use of this approach. We note that 
prev ious work on coalescing black hole-neutron star bina- 
ries (|Lee. Ramirez Ruiz fc Granoti[2005t lRosswoEjl2005|) 
has made use of t he p seudo-Newtonian potential of 
iPaczynski fc Wiital (|1980f ) for the black hole, which re- 
produces the existence and position of a last stable orbit 
for test particles in circular orbits around a Schwarszchild 
black hole. Doing this thus requires considering fairly 
small mass ratios, a condition which is not met in the 
present set of calculations, with the additional complica- 
tion of highly eccentric orbits for the parabolic encoun- 
ters of interest. We have hence elected to keep to a purely 
Newtonian formulation which, albeit simplistic, is quan- 
tifiably so. 

For neutron stars, most equations of state reveal that 
the radius varies little over a range of masses. In the 
case of polytropes, the mass-radius relation is R oc 
M (r-2)/(3r-4) ; g0 jf the ac ji a batic index is T = 2, the 
radius is only a function of the structure constant K and 
the central density, p c . Our standard "neutron star" is 
thus a spherical, non-spinning polytrope with T = 2, 
-Wns = 1-4M Q and i?NS = 13-4 km. In order to inves- 
tigate the effects of a different compressibility, we have 
also considered a neutron star with index V = 5/3 of 
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the same mass and radius for certain orbital parameters, 
detailed below. We have also considered the case of a 
low-mass white dwarf (Mwd — 0.5M Q ) interacting with 
a black hole. In this case it is appropriate to use the 
equation of state for a cold non-relativistic Fermi gas, 
P = ifn.r.p 5 / 3 , giving a radius i?wD = 1.1 X 10 4 km. 

Three-dimensional calculations of binary interactions 
are typically evolved for only a few tens of milliseconds, 
and during this short time large scale gravitational dy- 
namics determine the final state of the system. Once the 
initial conditions are set, we use a simple ideal gas equa- 
tion of state, where the pressure is given by P = pu(T— 1) 
and u is the specific internal energy, to follow the ther- 
modynamics of the gas, with no heating or cooling mech- 
anisms present. Shocks are allowed to form in the usual 
numerical way (through an artificial viscosity), and are 
the only source of local dissipation. 

The single additional ingredient that is necessary in 
terms of global dissipation is the emission of gravita- 
tional waves, since it can (and does) affect the orbital 
evolution. For binary mergers, the point-mass approxi- 
mation in the weak field limit is often used to compute 
an effective drag, removing angular momentum and en- 
ergy from the system. It is switched off once the stars 
come into contact, or if one of them is tidally disrupted. 
However, for the present set of calculations on various 
orbits it is not a good approximation. Rather, we re- 
quire an expression for the energy loss rate of extended 
bodies, which can be used for more general trajectories 
(for example initially parabolic). We thus compute the 
rate of energy loss as 



(ii?GW 1 G 



dt 



5 c 5 



dljk dljk 
dt dt 



(1) 



where Ijk is the traceless mass quadrupole moment. 
The first two time derivatives can be easily calculated 
without having to perform two numerical time deriva- 
tives expli citly with the use of the continuity and Eule r 
equations (jRasio fc Shapirolll992HLee fc Kluzniakll 19990 . 
The third derivative requires numerical treatment, and 
numerical noise in its calculation can be reduced by in- 
terpolation over neighboring points in the evolution. We 
have tested this implementation by computing the evolu- 
tion of binaries with large separations, where the orbital 
decay time scale is much longer than the orbital period, 
and find excellent agreement with the point-mass formula 
for circular orbits. It is important to consider this more 
accurate treatment for parabolic encounters, since the 
energy loss rate at periastron can be substantially larger 
(by a factor of 2-3, depending on the impact parame- 
ter) than that estimated from the expression for point 
masses in circular orbits, leading to different, and over- 
all, faster encounters. In addition, the effect on the loss 
rate because of the formation of tidal bulges on the stars 
is automatically taken into account, since it is the full 
fluid quadrupole moment that is used in computing the 
loss rate. 

2.2. Conditions for stellar encounters 

The stellar velocity dispersion in GCs is insignificant 
when compared to that acquired from the mutual gravi- 
tational acceleration of two compact, stellar mass objects 



as they approach each other. It is thus reasonable to as- 
sume at first that encounters will be parabolic, i.e., with 
vanishing orbital energy, E, at infinity. However the im- 
pact parameter, 6, may vary, and is related to the total 
orbital angular momentum, L, in the two-body system. 

When the separation between the two masses, a, is 
large compared to their individual radii, R ll we may 
safely consider that they behave as point particles and 
compute the orbit accordingly. It is only when these 
quantities become comparable that finite-size effects 
need to be considered, and a full hydrodynamic calcu- 
lation must be performed. It has been shown for the 
case of compact binary mergers that tidal effects can 
de-stabilize the orbit even in the absence of relativis- 
tic con si derations at small enou gh separations ()Lai et al.1 
Il993allbl ; iRasio & Shapiro 1994). The orbital decay rate 
thus induced is comparable to that due to gravitational 
radiation back reaction, and can even dominate the evo- 
lution for the last few cycles. The equation of state plays 
a role in the magnitude and growth of this effect, and for 
neutron stars it is almost certainly generically impor- 
tant. We have thus elected to study collisions by placing 
the stars (whether they are black holes, neutron stars 
or white dwarfs) initially at separations comparable to 
those used in merger calculations. This choice is always 
a trade-off, since using a larger separation will always be 
more accurate, at the cost of added computational effort 
(during which the components approach each other with 
little visible evolution). The point worth noting here is 
that as a result of tidal effects, a bulge will form on each 
star as they approach, but will slowly fall behind the line 
joining the centers. At periastron, a significant lag angle 
may have developed, and this can have an effect on the 
subsequent evolution of the material, through the action 
of gravitational torques by the more massive primary on 
the lighter, bar-like secondary. For the present set of cal- 
culations we have chosen to use spherical, non-perturbed 
stars for initial conditions. A solution in terms of com- 
pressible t ri-axial ellipsoids is possible for the previous 
evolution ([Lai et al.l 1199 3b). but given our overall sim- 
plifications we believe this will not affect our results sig- 
nificantly. At the chosen initial separations, the induced 
tidal bulges and lag angles are still quite smal0. 

A convenient way to parameterize the encounters is by 
comparing the strength of the tidal field produced by the 
primary to the self gravity of the secondary, through 



V 



Mi R\ 



1/2 



(2) 



where M\ and M 2 are the primary and secondary mass, 
R p is the periastron distance, and R2 is the radius of 
the secondary. With this definition, at a fixed mass ra- 
tio q = M2/M1 < 1, encounters with rj ^> 1 have large 
impact parameters, while rj ~ 1 marks collisions where 
the secondary is tidally disrupted because the tidal field 
has become as intense as that holding the star together. 
An important difference between encounters with com- 
ponents of similar mass and those in which q <C 1, as 

2 A detailed comparison of the differences between the use of 
spherical stars as oppose d to tri-axia l ellipsoids was carried out in 
the case of mergers (Lee 2000, 2001). Overall, results differed by 
less than RJ 5% for the equation of state used here with T = 2. 
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for the disruption o f stars by supermassive black holes in 
AGN (iFrankJ U97l ILacv et al.lll98l ICarter fc LuniTnetl 
119831 : lReedll988f> . is that the periastron distance is much 
smaller in the former, and in fact comparable to the stel- 
lar radius. For a given intensity of the encounter, the 
periastron distance can be written as 

R P =R 2 V V3 q- 1/3 - (3) 

Thus for a disruptive encounter with r\ ~ 1, R p ~ 200i?2 
if q = 1CT 7 , but R p ~ 2R 2 for q ~ 1/3. In the first 
case the disruption is a purely gravitational encounter, 
while in the second it is additionally a direct collision, 
which may modify the mass ratio substantially during 
the encounter. 

The parameter rj can also be thought of as the ra- 
tio between the dynamical time scale of the star (or its 
rotational break up period) and the duration of the en- 
counter. Thus for the intrinsic spin of the secondary to 
be of any relevance, it must be rotating near the Ke- 
plerian limit. The most rapidly spi nning neutron stars 
detected in LMXBs in our gala xy (|Backer et al.l 11982ft 
and in GCs (jHessels et al.l I2006h are clearly below this 
threshold, by a factor 3-5, depending on the assumed 
equation of state. It is thus reasonable to assume that 
in most cases the neutron star spin is negligible for the 
purposes of the encounters considered here, and we do so 
in what follows. 

2.3. Tidal capture and the formation of close binaries 

When two compact stars pass close to each other so 
that ?? — 1, the tidal perturbations produce accelerations 
of the stellar material relative to the stellar center. The 
internal energy gained by the stars is taken from the 
kinetic energy of their relative motion. While the internal 
energy AEt gained by each component is a small fraction 
of the star's internal energy, it can be comparable with 
the kinetic energy of relative motion of the two bodies at 
large separation. The two stars will thus become bound 
if enough energy is absorbed in exciting the oscillations. 
The resulting elliptical orbit immediatey after capture 
will have an initial eccentricity only slightly less than 
unity. However, there will be many successive passages 
at about the same periastron distance, and these will lead 
to further energy transfer. Ultimately the orbit should 
become nearly circular. 

The value of AEt and thus the resultant condition for 
tidal capture may be determined by requiring that the 
frequency of the perturbation be slightly smaller than the 
natural frequencies of the perturbed system. The tidal 
field excites non-radial modes of oscillation in the sec- 
ondary and, while there are many of these, those most 
effectively excited have the smallest number of nodes, 
with an angular frequency uj osc that is usually some two 
to three times oj 2 = {GM 2 / R^) 1 ^ 2 ■ Now, the angu- 
lar frequency of the encounter u> p is effectively V p /R p 
and can be written as lu p w 2lu 2 (R 2 / i? p ) 3 / 2 . Thus for 
R P /R 2 — 2, uj p /u osc w 0.3 for the lowest modes. How- 
ever, for R p /R 2 = 10, lu p /uj osc w 0.03. For such slow 
changes, the shape of each star will adjust to the chang- 
ing form of the equipotential surfaces and the net heating 
is markedly reduced. The binaries formed by tidal cap- 
ture are thus generally very hard. 

Since the fraction of the initial angular momentum 
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Fig. 1. — The encounter of a neutron star with a black hole at 
relatively large impact parameter (run L3) leads to a large number 
of periastron passages before tidal disruption. Here we show the 
separation (in km) between the neutron star and the black hole 
from the start of the simulation until full disruption of the star, 
taking about one quarter of a second. 

which is transferred to stellar rotation can scarcely ex- 
ceed a few percent because of the relatively small stellar 
radius of gyration, one can assume that the orbital an- 
gular momentum remains roughly constant. Hence in 
a circular orbit of radius R c and relative velocity V c , 
R C V C must equal the initial R p V p . With V p obtained 
by the condition that (l/2)^ r U p 2 = GM 2 2 /i? p , where 
/!,. = M\M 2 j [M\ + M 2 ) and V c computed from force 
balance in a circular orbit, one finds R c = 2R p . To ob- 
tain more exact results for the resulting conditions for 
tidal capture, the excitation of the individual normal 
modes must be considered, and the h eating summed over 
all mo des. Analysis of this effect by iPress fe Teukolskvl 
(|1977f ) has taken the parabolic motion of the two stars 
accurately into account, but in the tidal potential only 
terms varying as 1/r 2 and 1/r 3 have been considered, 
an approximation valid for R p /R 2 > 3. Detailed cal- 
culations were first carried out using linear theory for 
r = 4/3 (IPress fe Teukolskvlll977t ILee fe Ostrikenll98l 
and for r = 5/3 (|Lee fe Ostrikeri Il986f) . They were 
follo wed by many other studies using bot h linear the- 
ory ([McMillan et al.lll987t lKochanekjll992ft and numer i- 
cal hydrodynamical calculations (jRasio fc Shapiroll991| ) . 
Such treatments give R p /R 2 between 2 and 3 as a con- 
dition for tidal capture (the exact value depending on 
the properties of the stellar structure model used). Thus 
an approximate condition for tidal capture may be given 
as R p < 3i?2 irrespective of the precise details since for 
R P /R 2 < 3, uj p ~ uj osc . 

For compact binaries, angular momentum losses to 
gravitational waves, as well as dynamical mass ejection 
from the system must additionally be considered. The 
results for a black hole-neutron star encounter with mass 
ratio q = 0.31, 77 = 3 (corresponding to R p /R 2 = 3.1), 
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TABLE 1 

Parameters for orbital encounters. 



Run 


Prim. 


Sec. 


r 


Mi 

M 


M 2 

M 


M 2 
1= Ml 


V 


Hp 
«2 


Lo 


BH 


NS 


2 


4.51 


1.4 


0.31 


1.0 


1.5 


L or 5 / 3 


BH 


NS 


5/3 


4.51 


1.4 


0.31 


1.0 


1.5 


Li 


BH 


NS 


2 


4.51 


1.4 


0.31 


1.5 


1.9 


L 2 


BH 


NS 


2 


4.51 


1.4 


0.31 


2.0 


2.3 


L 2 r 5/ 3 


BH 


NS 


5/3 


4.51 


1.4 


0.31 


2.0 


2.3 


L 3 


BH 


NS 


2 


4.51 


1.4 


0.31 


3.0 


3.1 


NSNS 


NS 


NS 


2 


1.75 


1.4 


0.80 


1.0 


1.0 


BHWD 


BH 


WD 


5/3 


2.5 


0.5 


0.20 


1.0 


1.7 



where the neutron star is modeled as a polytrope with 
r = 2 (run L3 in Table [TJ are shown in Figure [TJ As 
expected, the strong tides trigger complex oscillations 
in the secondary star which can be clearly seen in the 
variation of the maximum density in the core after the 
first periastron passage. Further energy transfer occurs 
in subsequent passages as the orbit becomes tighter and 
progressively circular, until eventually some direct mass 
transfer takes place along with mass stripping. We find 
that the secondary is not shredded immediately, but only 
after more than a dozen passages at about the same pe- 
riastron distance. Our condition for for tidal capture is 
in rough agreement with i? p /i?2 < 3 although the in- 
clusion of gravitational waves as well as matter ejection 
from the system allows for a more rapid variation of an- 
gular momentum disposal. For 77 ~ 1, the secondary star 
will experience a direct physical collision at least of its 
outer layers. The process is rather complex with some of 
the gas escaping entirely from the system on outbound 
trajectories. The core may preserve its integrity for a few 
orbital periods around the primary in a common enve- 
lope before finally coalescing. It is to this problem that 
we now turn our attention. 

A summary of the most important aspects of the ini- 
tial conditions thus chosen is given in Table [TJ There 
are six parabolic encounters involving neutron stars with 
black holes, one double neutron star collision, and one 
black hole- white dwarf encounter. The parameter 77 has 
been computed from the point mass orbital parameters 
in Newtonian gravity. Naturally the centers of mass of 
each star do not follow these solutions because of finite 
size effects and the emission of gravitational waves, but 
they allow for a characterization of each case. Note that 

since r\ oc Rp^ 2 and L oc R^/ 2 , the change in 77 by a factor 



3 is actually a variation of 3 1 / 3 
momentum. 



1.4 in orbital angular 



2.4. Stellar disruption and disk formation 

For binary mergers, tidal disruption usually occurs af- 
ter a single periastron passage, because the system is 
already very tightly bound by gravity. In the case of 
collisions, we find that the secondary is not shredded im- 
mediately after the first interaction. It does, however, 
lose a substantial amount of orbital energy and angu- 
lar momentum through two main channels: emission of 
gravitational waves and transfer by gravitational torques 
through the formation of the tidal bulge. The bulge in 
effect deforms the spherical star into a bar, which the pri- 
mary can then spin-up during the brief encounter. The 
strong tides additionally trigger complex, nearly radial 




Fig. 2. — The neutron star core mass decreases during the colli- 
sion with a black hole by tidal mass stripping. As the strength of 
the encounter, r\, decreases, the neutron star survives for a greater 
number of periastron passages before being eventually shredded. 
The secondary's mass (in solar masses) is shown for runs Lo, 
LoT5/3, Li, L2 and L2r 5 /3 until disruption. The thick (thin) lines 
are for V = 2 (r = 5/3). 

oscillations in the secondary which can be clearly seen in 
the variation of the maximum density in the core after 
periastron passage. Some direct mass transfer occurs as 
well, along with mass stripping which can form a first 
accretion structure around the primary, which is more 
massive for low initial impact parameters. Additionally, 
some matter is flung out through the exterior Lagrange 
point to large distances, creating a tidal tail. 

In runs Lo and Li, the neutron star core (which now 
contains ~ 1 M Q and 1.2 M respectively, see Figure [5]) 
does not survive the second encounter, and forms a mas- 
sive disk around the black hole, as well as an elongated 
tidal tail. In run L2 there is enough energy and an- 
gular momentum at the outset that the core is able to 
return a third time, after which it too is fully disrupted 
and forms a disk (see Figure [3]) . Each successive passage 
feeds the accretion disk and simultaneously forms a tidal 
tail, which is not set to collide with previous ejections. 
We show in Figure HJ the time evolution for different ini- 
tial impact parameters, projected onto the orbital plane. 
Even with three periastron passages, the collision and 
disruption are essentially over after ~ 50 ms, because 
the initial passage drains enough orbital energy and an- 
gular momentum to bind the system very effectively. 

For runs carried out with a soft equation of state, 
r = 5/3, the results are qualitatively the same for low 
impact parameter, with the main difference being in the 
spatial extent of the accretion disk formed, and in that 
of the tidal tails, both being greater than for r = 2. For 
large impact parameter, the gravitational interaction is 
quantitatively different during the first periastron pas- 
sage, since the star effectively resembles a point mass to 
a greater degree in run L2r 5 / 3 . Angular momentum and 
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Fig. 3. — The separation (in km) between the neutron star and 
the black hole is shown until disruption of the neutron star for runs 
Lo, Lor 5 /3, Li, L2 and 1^2^/3. Note the change in scale on the 
time axis when compared with Figure Q] The thick (thin) lines are 
for T = 2 (r = 5/3). 
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Fig. 4. — Parabolic collisions of neutron stars with black holes. 
Logarithmic contours of density in the orbital plane (equally spaced 
every 0.25 dex) are shown for runs Lo (left), Li (middle) and L2 
(right). The lowest contour in bold is at logp [g cm -3 ] = 10. The 
time elapsed since the beginning of the simlulation is indicated 
in each panel (increasing top to bottom), as is the distance scale. 
Note the different number of periastron passages and tidal tails 
formed until final disruption in each case. 
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Fig. 5. — The total angular momentum as a function of time is 
shown for runs Lo, Lor 5 /3, Li, L2 and L2r 5 /3, in units of L* = 

Mtot VG^Rns 1 where Mtot is the total mass and fi is the reduced 
mass. For each case, successive periastron passages are marked 
on the curve as Pi, P2, ... .The thick (thin) lines are for T = 2 
(T = 5/3). 

energy transfer through torques is thus less efficient and 
the core of the neutron star is transferred to a higher 
orbit than in run L2 (note the difference in secondary 
apocenter values in Figure [3]). At the secondary passage 
l__ 26 ms after the start of the simulation, however, the 
star is fully shredded by tidal forces and the final disk 
forms promptly. 

For the disruption of a star by a supermassive black 
hole, the former moves essentially in the fixed back- 
ground metric imposed by the hole, and this allows for 
a simplified treat ment of the dynamics. It was found 
in earlier studies (Rccs 1988|), both through analytical 
considerations and direct numerical simulation, that es- 
sentially half the mass of the star is dynamically ejected, 
while the remaining half is captured by the black hole, 
being on eccentric trajectories which, with variable delay, 
will bring them back to the vicinity of the primary. 

There are two important differences between such a 
scenario and that considered in the present study. First, 
as already noted, the mass ratio is of of order 0.1-1 in- 
stead of 10 -7 , which can produce a direct collision with 
accompanying mass transfer simply by virtue of the small 
periastron distance of a disruptive event. Second, at the 
small distance scales implied by the fact that we are con- 
sidering compact stellar mass objects, gravitational wave 
emission is intense, and can drain a substantial fraction 
of the total kinetic energy and angular momentum dur- 
ing a single passage (see Figure [SJ. Thus the system 
becomes non-conservative from a point of view of orbital 
dynamics, and the previous reasoning does not apply. 

After the neutron star has been fully disrupted, the 
remnant consists of a black hole surrounded by a torus, 
and a series of tidal tails, depending on the number of pe- 
riastron passages which occurred. The tori are typically 
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Fig. 6. — Logarithmic contours of density in the orbital plane 
(equally spaced every 0.25 dex) for the collision of two neutron stars 
(run NSNS) on a parabolic orbit with initial mass ratio q = 0.8. 
The lowest contour in bold is at logp [gem -3 ] = 10. After an 
initial passage the secondary is entirely disrupted and wraps around 
the primary, forming an envelope and a single tidal tail. 

200-300 km across and contain M^ is ^ ~ O.IMq. They are 
comparable in size and mass to those encountered during 
binary mergers, with densities p ~ 10 11 — 10 12 g cm -3 
and internal energies u ~ 10 18 — 10 19 erg g _1 , equivalent 
to 1-10 MeV/baryon. The binary interaction is violent 
and complex, though, and this can be qualitatively seen 
in the fact that by the end of our calculations (roughly 
20 ms) the accretion structures still show significant de- 
viations from azimuthal symmetry. 

For double neutron star encounters, the mass ratio is 
even greater (we have computed one collision for q = 0.8, 
which is likely to be a lower bound for such systems). For 
the case considered with adiabatic index T = 2, the two 
stars are actually the same size, and directly impact each 
other for an encounter strength rj = 1. The initial colli- 
sion binds the stars in an elliptical orbit but does not lead 
to the formation of a significant tidal tail. The less mas- 
sive secondary is strongly distorted and spun up, and a 
bridge of material temporarily joins the two stars. Upon 
a second passage the secondary is entirely shredded and 
wraps around the primary (see Figure [BJ . The material 
from each stars remains largely separate, with that from 
the primary remaining essentially in the core, in a man- 
ner similar to what is obtained for mergers of unequal 
mass neutron stars. 

The final remnant in this case consists of a slightly dif- 
ferentially rotating core of 2.4M Q with radius i? C orc — 
20 km and a maximal rotation frequency of f2 max — 
3000 rad/s, corresponding to a period of 2.1 ms (see 
Figure [7]). This is surrounded by an envelope out to 
~ 80 km containing M cnv ~ O.49M0, where the orbital 
frequency is Q oc r . In the outer tail of material the 
rotation is essentially Keplerian, with £1 oc r~ 1,5 (the ro- 
tation profile is shown as well in Figure [7]) . Note that 
the core is above the threshold for collapse of a cold, 




r Tkml 

Fig. 7. — Azimuthally averaged angular velocity profile (solid 
line) and enclosed mass (dashed line) in the double neutron star 
collision remnant. The rapidly rotating inner core (P=2.1 ms) is 
surrounded by a large envelope. The rotation of the tidal tail is 
visible at radii greater than 100 km, and reference power laws are 
given for the envelope and tail. 

non-rotating configuration in most equations of state, 
but could avoi d this given the rapid and differential ro - 
tation present (|Cook et al.lll994HBaumgarte et al.ll2000n . 
In addition significant heating of the core can raise the 
critical threshold mass for collapse to ~ 1.35M co id 
(Shibata & Taniguchi 2006). In this particular case the 
core could conceivably remain stable for a longer time, 
spinning down on a secular time scale due to the emis- 
sion of gravitational waves or magnetic torques. If it were 
to indeed form a black hole, the distribution of angular 
momentum in the envelope is such that only material be- 
tween 50 and 100 km, amounting to 0.17M Q , would have 
enough centrifugal support to form an accretion disk. 

For completeness we have also considered the interac- 
tion of a low mass (Mwd = 0.5M Q ) white dwarf with a 
compact object (M co = 2M@). The latter could be either 
a massive neutron star or a low mass black hole. Numer- 
ically it is irrelevant because the characteristic scales of 
the two objects are so different that it is impossible to 
resolve simultaneously the black hole horizon (or neu- 
tron star surface) and the entire white dwarf. For the 
actual calculation the absorbing accretion boundary has 
been placed 100 gravitational radii from the center of the 
compact object. 

The white dwarf is modeled as a polytrope with 
adiabatic index T — 5/3, appropriate for a cold, 
non-relativistic degenerate configuration, and the pres- 
sure is thus given by P — K n . r .p 5 / 3 , where K n , T , — 
(3/7r) 2 / 3 ft, 2 /(20TO e (2m p ) 5 / 3 ). With the given mass, the 
stellar radius is i?vvD = 1.1 X 10 9 cm, a factor / ~ 800 
larger than our standard neutron star. In fact with the 
same input physics one could simply scale the results 
from the BH-NS interaction (given an identical mass ra- 
tio) by increasing distances by a factor / and temporal 
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Fig. 8. — Logarithmic density contours (equally spaced every 0.25 
dex) for the encounter between a low mass white dwarf and a com- 
pact object. The lowest contour in bold is at logp [g cm -3 ] = 3. 
Note the difference in scales (spatial and temporal) when compared 
to those in Figures [4] and [6] The higher compressibility of the ma- 
terial also makes for a wider distribution of ejected material and a 
generally wider distribution of gas. 



TABLE 2 

Disk, tail and ejected masses. 



Run 


M diBk /M e 


M ts , a /M Q 


Mej./iW© 


Lo 


0.23 


0.15 


0.05 


L r 5/3 


0.25 


0.11 


0.04 


Li 


0.16 


0.31 


0.13 


L 2 


0.09 


0.47 


0.19 


L 2r 5 /3 


0.11 


0.19 


0.04 


NSNS 


0.17 


0.13 


0.03 


BHWD 


0.13 


0.35 


0.24 



scales by / 3 / 2 ~ 2.3 x 10 4 ( giving hundreds of seconds 
instead of tens of milliseconds). This is strictly not cor- 
rect, however, since gravitational radiation reaction in- 
troduces an absolute scale into the problem, and energy 
and angular momentum losses through this channel are 
insignificant in comparison to the BH-NS case. The in- 
teraction proceeds then at a more leisurely pace, so much 
so that we were unable to follow it to a second perias- 
tron passage, even though the simulation covered nearly 
one thousand seconds. At this stage stripping during the 
close passage has formed a torus around the compact ob- 
ject, linked to the stellar core by a long and narrow bridge 
of material (see Figure [5]). The core itself is rapidly spin- 
ning due to gravitational torques exerted by the primary, 
and the typical tidal tail has formed at large radii. 

Since we are considering encounters in which disrup- 
tion occurs practically by construction (given the choice 
of the parameter 77), the gravitational torques exerted 
on the secondary are of comparable magnitude in all 
cases (with appropriate scalings). We find that the an- 
gular frequency of the core after the first periastron pas- 
sage is ~ (0.2 - 0.4)fi , where fl = [GM 2 /Rl) 1/2 is 
the break up rotation frequency of the unperturbed sec- 
ondary. Thus the spin periods are ~ 2 ms and ~ 8 s 
for encounters involving neutron stars and white dwarfs 
respectively. 

The mass of the disk (see Table [2]) present at the 
end of the calculation as a result of the encounter is 
computed from the fluid elements in close, essentially 
circular orbit about the central object. It is typically 
Mdi s k — (0.1 — 0.3)M Q , regardless of the type of en- 
counter. Likewise, the mass of the tidal tails, M ta ii, 



Fig. 9. — Tidal tails formed by the disruption of a neutron star 
by a black hole (run Li) and a white dwarf (run WDBH). For 
the neutron star disruption the impact parameter was such that 
two periastron passages occurred, giving rise to two distinct ejec- 
tion structures thousands of kilometers across by the end of the 
simulation. For the case involving the white dwarf, only one pas- 
sage has occurred and the core of the star is still present. Note the 
much larger spread in the distribution of the fluid, due to the lower 
adiabatic index used in the case of the white dwarf. 

is obtained by adding over all the fluid elements within 
these structures, whether they are bound to the central 
mass or not. For black hole- neutron star encounters there 
is a clear trend of decreasing disk mass and increasing 
tail mass as the initial impact parameter grows at fixed 
compressibility. The total amount of matter dynamically 
ejected (computed as that with positive total energy) also 
increases, and is a result of the greater number of ejection 
episodes associated with larger values of the orbital an- 
gular momentum. For run L2 we find that nearly 0.2M@ 
are thus lost. For the double neutron star collision the 
characteristics of the system (total mass and mass ratio) 
are such that the ejected mass is significantly lower, al- 
though a substantial tail is still present. We will return 
to the implications of mass ejection below. Including the 
effects of General Relativity will likely alter these values, 
by up to one order of magnitude if the differences be- 
tween calculations performed in Newtonian theory and 
those using GR for mergers are used as guidance. 

2.5. Tidal tails and mass ejection 

For the two members of the system to come together 
and eventually merge, they must lose energy and angular 
momentum. This can be achieved through the emission 
of gravitational waves, or the ejection of matter, or a 
combination of both. In either a binary system in circu- 
lar orbit or a parabolic approach there is a substantial 
amount of rotation. Since a small amount of matter, re- 
moved to a large radius, can carry a great deal of angular 
momentum, the formation of tidal tails out of material 
stripped through the outer Lagrange points in the effec- 
tive potential is an efficient way to produce a single object 
at the center. This is the fundamental reason why such 
structures form in either stellar or galactic collisions. 

Just as the vanishing orbital energy of a parabolic or- 
bit implies successive periastron passages as the core be- 
comes more and more bound to the primary, so too it 
is easier to dynamically unbind matter to infinity when 
compared to a bound binary. As each passage proceeds a 
new ejection episode occurs, giving rise to an additional 
tail (Figure [5]) . The amount of mass thus ejected, and 
shown in Table [5J is considerably greater than for a bi- 
nary coalescence, by about one order of magnitude (|Led 
[20011) . 

In all cases when multiple ejection events occur (we 
observed up to three for run L2), the velocities and orien- 
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tations are such that the initial tail will not be overtaken 
by latter ones. This is simply because the first one has 
a clear head start, but also because subsequent events 
occur at different orbital phases. 

An interesting point regarding the tails concerns their 
hydrodynamic stability. Their motion is essentially bal- 
listic, dominated by the potential well of the central 
mass. They are nevertheless susceptible to the vari- 
cose, or sausage instability firs t identified by Rayleigh 
in f899 (|Chandrasekharl 1196 If ). This is due to self- 
gravity and occurs for cylindrical configurations of an 
incompressible fluid for perturbations with wavelength 
A > A* = 2irR cy \/x*, where x* ~ 1 and i? cy i is the 
radius of the cylinder. The fastest growing mode has 
x = 0.58, wavelength A ~ lli? cy i and a growth time 
T = A/iA-KGp) 1 / 2 (|Chandrasekharl[l96l . For sufficiently 
stiff equations of state (with adiabatic index T > 2.5), 
even though not strictly incomp ressible, this is actually 
seen in n umerical simulations (Ra sio &: Shapiro] Il994t 
iLeei [20 00): on a time scale given approximately by r, 
condensations form, and are separated roughly by the 
wavelength of the fastest growing mode given above. For 
softer equations of state (even such as the one used here 
for neutron stars when r = 2) the effect is not present, 
and even less so for the calculation involving the neutron 
star with T = 5/3 and the white dwarf (Figure [S]). 



3. EMISSION OF GRAVITATIONAL WAVES 

In the case of merging binaries, the early gravitational 
waveforms, when the separation is much larger than the 
stellar radius, can be computed using the weak field 
approximation analytically, and from them the stellar 
masses may be accurately determined. As the stars be- 
come distorted by the tidal field, the signal deviates form 
this solution and finite size effects accelerate the decay. 
The secondary is then fully accreted by the black hole 
(if one is present) or tidally disrupted, and the emission 
abruptly ceases. The precise frequency where this occurs 
can lead in principle to accurate deter minations of the 
neutron star radius (jFaber et al.l 120021) . and thus, since 
the mass is already known, to useful constraints on the 
equation of state at supra-nuclear densities. 

The total mass and relative velocities involved in a 
parabolic encounter are similar to those encountered in 
close binaries, and so we would expect the characteris- 
tic frequencies and strength, or amplitude of the signal 
in gravitational waves, to be comparable in this case. 
The biggest difference arises, and is crucial in terms 
of detectability, because the collision does not involve 
a leisurely spiral in over many orbital periods, and as 
thus lacks the characteristic precursor chirp signal which 
would slowly sweep through the interferometer's detec- 
tion band. Figure [10] shows the computed waveforms 
(one polarization is given) for black hole-neutron star 
encounters and the double neutron star collision. 

A neutron star (M2) approaching a stellar mass black 
hole (Mi) with r\ ~ 1 will be disrupted in a single pas- 
sage and the particles in the disrupted remnant follow 
approximately independent Keplerian orbits. The de- 
tectable gravitational signal will thus have a burst-like 






Fig. 10. — Gravitational waves emitted as a function of time for 
runs Lo, Lor 5/ /3, Li, L2, L2r 5 /3, and NSNS as seen by an observer 
placed along the rotation axis (r/14. is plotted (in m), where r is the 
distance from the source to the observer). The thick (thin) lines are 
for r = 2 (r = 5/3). The vertical range is identical in all frames, 
but note the different scaling on the time axis. In all runs except for 
Lo, Lor 5 /3, and L2r 5 /3, the oscillations of the neutron star core 
following the first and even second periastron passage are clearly 
visible. When varying the stiffness of the equation of state, the 
waveforms are nearly indistinguishable for small impact parameter 
and clearly separated for large values, due to the difference in the 
associated orbital evolution. 

behaviour, roughly characterized by an amplitude 
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Here M denotes the total mass, Mi + M 2 . LIGO will 
be able to detect gravitational wave from impact involv- 
ing neutron stars and stellar mass black holes if r] < 1 
and the distance is D < 50 Mpc. In what follows, we 
compare these simple estimates of gravitational radiation 
against more detailed results obtained with our numeri- 
cal scheme. 

The signal exhibits a local peak at each periastron pas- 
sage, until the time when the star is completely disrupted 
and the amplitude vanishes (note the different time scales 
on each plot in Figure [TU)) . In runs Li, L 2 and NSNS, 
smaller oscillations of decaying amplitude after each pas- 
sage (but perhaps the last) are also clearly present when 
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Fig. 11. — Gravitational wave energy spectrum for runs Lo, 
Lor 5 /3, Li, L2, L2r 5 /3 and NSNS. The ranges are identical in 

all frames. The reference power law is dE/df oc f -1 ^ 3 , the char- 
acteristic spectrum for the in-spiral of a point mass binary. The 
thick (thin) lines are for T = 2 (r = 5/3). The oscillating nature 
of these spectra is due to the finite extent of wave trains present 
in the amplitudes during different segments of temporal evolution. 

r = 2. For instance, for run L2 between 3 and 15 ms, 
their frequency is ^ osc ~ 1750 Hz. This is not due to the 
rotation of the neutron star (the spin frequency of the 
core is at this point approximately i^pin — 320 Hz) but 
essentially to radial vibration modes excited by the ac- 
tion of the primary at periastron passage. The frequency 
of such modes is close to the natural value given by hy- 
drostatic equilibrium, namely v osc — (GM/ i? 3 ) 1 / 2 /2tt ~ 
1400 Hz. The small impact parameter and rapid disrup- 
tion of the star in run Lo, and the high compressibility 
used in run L2r 5 / 3 does not allow for a clear manifesta- 
tion of such oscillations. 

The power spectrum of these gravitational waves is 
shown in Figure [TT] There are peaks at 1900, 1600 and 
1000 Hz for runs Lo, Li and L2, respectively. These 
correspond to the interaction time for each run, which 
increases (thus decreasing the frequency) as the impact 
parameter (or equivalently, the total angular momentum 
in the system) increases. Secondary periastron passages 
make this peak wider at higher frequencies. This is par- 
ticularly evident when comparing the results of runs L2 
and L2r 5 / 3 . The latter has a single late-time perias- 
tron passage (instead of three for the former), and the 
peak is broadened to ~ 2000 — 3000 Hz. We note that 
the oscillation frequency of the core, visible in the sig- 
nal when r — 2, is undetectable in these spectra, being 
buried within the main peak due to the binary inter- 
action (this is to be expected, since by definition of a 
tidally disruptive event, the interaction time is roughly 
equal to the free fall time scale of the star). The spec- 
trum for run NSNS shows a similar overall structure. 
We note that it is difficult to extract a clean spectral 
signal with a high signal to noise ratio, because the sim- 




Fig. 12. — Power radiated in gravitational waves for black hole- 
neutron star encounters and the double neutron star collision. 
There is one successively weaker peak for each periastron passage 
until complete tidal disruption. The curve labeled "Binary" shows 
the luminosity computed for a black hole-neutron star binary with 
the same mass ratio, q = 0.31 and initial separation, ri = 3.7-RnSi 
as runs Lo through L2 , for whic h disruptio n occurs promptly after 
the beginning of mass transfer (Lec 2001). The thick (thin) lines 
are for T = 2 (r = 5/3). When varying the stiffness of the equa- 
tion of state, the luminosities are nearly indistinguishable for small 
impact parameter and clearly separated for large values, due to the 
difference in the associated orbital evolution. 

ulated time is quite limited. Secondary high frequency 
variations at ~ 5000 Hz can be seen in the spectra for 
runs Lo, Li and NSNS and are related to the asym- 
metry of the waveform around the time of closest ap- 
proach. The large number of periastron passages in 
run L2 makes the spectrum more noisy and this fea- 
ture is not observed. Finally, we also show the grav- 
itational wave luminosity in Figure 1121 the integration 
of which reveals that the efficiency for their emission is 
e = E GW /(Mc 2 ) = 10~ 2 , 7 x 10~ 3 , 6.8 x 10" 3 , 5.9 x 10~ 3 
for runs Lo, Li, L2 and NSNS respectively. For compari- 
son, the efficiency in the case of a merging binary with the 
same mass ratio is « 5 x 10 -3 (|Ledl206 l). We note also 
from Figure rT2] that the maximum power radiated during 
a collision and even during subsequent close approaches 
can be significantly higher than that in a merging binary 
with the same mass ratio, due to the larger velocities in 
an eccentric orbit at periastron. 

4. THE RATES OF COLLISIONS AND TIDAL CAPTURES 
OF COMPACT OBJECTS IN GLOBULAR CLUSTERS 

The highest stellar densities in GCs are reached in the 
core, in particular during core collapse. In this section 
we estimate the collision and tidal capture rate as a func- 
tion of time by first reconstructing the core evolution of 
a typical post-core-collapse GC. We then re-scale the re- 
sulting encounter rate evolution with the measured mass 
and size distribution of GCs in galaxies to obtain the 
rate per host galaxy. Finally, for a given galaxy luminos- 
ity density distribution, we obtain per volume of space 
and as function of time, the expected global rate of colli- 
sions and tidal capture of compact objects in the cores of 
GCs. The details of the encounter rate calculation will 
be presented in an forthcoming paper (van de Ven et al., 
in prep.), of which we give here a summary focussing on 
close encounters between compact stars. 

4.1. Close encounters 
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TABLE 3 
Relative encounter rates. 



type 1 type 2 M 2 /M R 2 /km £ col ^tid/^col -R min /i? 2 

NS NS 1.4 10 1.00 6.68 13.37 

BH BH 4.5 13 4.18 6.80 13.60 

WD WD 0.5 11000 392.87 2.99 5.99 

BH NS - - 4.85 7.18 16.51 

NS WD - - 747.11 7.12 7.12 

BH WD - - 1966.61 10.03 10.04 



rate follows from multiplying the (default) NS-NS colli- 
sion rate by the factor £ co i in column 5. In turn, multiply- 
ing the resulting collision rate with the factor in column 
'6 yields an estimate of the tidal capture rate. Here we 
■neglect the weak dependence of i? m in in equation (JSJ on 
tr c , which in turn varies only mildly during the evolution 
of the core. While in general a tidal capture does not 
neces sarily lead to the coal escence of the two stellar ob- 
jects ()Lee fc Ostriker| [l986). in case of compact objects 
.the approach is so close (see i? m i n / i?2 in column 7) that 
they will merge well within a Hubble time. In particular, 
for two neutron stars the close encounter rate is boosted 
by a factor ~ 6.7 due to tidal capture, without any sig- 
nificant delay with respect to collisions because of the 
very efficient emission of gravitational waves. 

The encounter rate might be significantly enhanced 
due to interactions between single stars and stellar bi- 
naries (and even between binaries), which have a much 
larger cross section, with i? m i n proportional to the bi- 
nary separation. Moreover, due to mass segregation the 
more massive stellar objects, including binaries, sink to- 
wards the center and can be important or even dominant 
in the core. Still, only the stars in very "hard" binaries 
with small enough separations (and thus smaller cross 
sections) might coalesce within a Hubble time. We dis- 
cuss binaries in Section 03] below, in particular focussing 
on compact binaries as they are commonly believed to be 
the progenitors of SGRBs. We restrict however the en- 
counter rate calculations to collisions and tidal captures 
between single stars, and hence consider the resulting 
values as lower limits. 

4.2. Evolving encouter rate 

After gradual contraction during an early phase, the 
core of a GC can go into d eep self-similar collapse 
(|Lvnden-Bell fc Eggletonl 119801) . The collapse is halted 
due to engery release from interaction s with binaries, also 
known as binary burning (Hills 1975), and/or from other 
sources includi ng a possible intermediate-mass black hole 
([Shapiro! 119771 ). The post- collapse evolution has bee n 
extensively studied (see e.g. Hcggic fc Ramamanil fl989'). 
but many aspects are still intensivel y investigated and 
debated. For example, iFregeaul (|2008|) recently proposed 
that the observed post-core-collapse GCs (~ 20% in 
the Milky Way) are most likely in the binary burning 
phase, while the remaining GCs are still in the process 
of core contraction and have not yet reached the binary 
burning phase. In the meantime, mat ching Monte Carlo 
models of M4 (IHeggie fc Giersd 120081) and of NGC 6397 
(|Giersz fc Hcggic 2009) seem to reveal that both GCs are 
past collapse in the binary burning phase, even though 
only NGC 6397 is one of the observed post-core-collapse 
GCs, while the surface brig htness of M4 is well fi tted 
by a standard King profile. IGiersz fc Heggid (|2009|) ar- 
gue that the differences in the surface brightness pro- 
files are most likely due to fluctuations in the core af- 
ter collapse. They point out that it is well possible that 
most GCs are post-core-collapse, but only a fraction hap- 
pens to be at the "peak" of the fluctuations that corre- 
sponds to a cusped surface brightness profile and leads 
to post-core-collapse classification. These fluctuations 
are not necessarily the well-known g r avothermal osci l- 
lations (|Bettwiese'r fc Sugimotol Il984t iGoodmanl I1987D , 
but could be the result of the stochastic nature of binary 



We assume the different types of stellar objects, i, are 
distributed homogeneously within a spherical core of ra- 
dius r c , each with fractional number /i < 1 and total 
number density n c . We further assume that the stars 
follow a Maxwellian velocity distribution function with 
dispersion a c . Together with the dominating gravita- 
tiona l focusing, this means we can approximate (Hcggic 
119751 ) the total collision rate as 

v co l = 2.1 x HT 3 Gyr- 1 -fi^x 

1 + 0i2 

V10 6 pc- 3 / VO.lpc/ ViOkms" 1 / X 
/ Mi + M 2 \ ( R min \ 
In^AlUtan")' (6) 

where <5i 2 — 1 if type 1 and 2 are equal, and 5\2 = oth- 
erwise. The stars have masses Mi, and their separation 
at closest approach is given by i? m i n , which we take to 
be the sum of the stellar radii, i.e., i? m in = Ri+Ri- 

In addition, we also consider encounters in which stars 
pass close enough to each other to form a binary by 
transfering orbital energy to internal stellar oscillations. 
We adopt a cross sec tion for tidal capture of the form 
(ILee fc Ostrikerl[l986l) 

where v- m i is the relative velocity at infinity, and u +i 2 = 
(2GAf 2 /i?2) 1 ^ 2 is the escape velocity at the surface of the 
se condary, captured star. We use the fitting functions 
of iKim fc Led (|1999l ) to obtain the amplitude a for en- 
counters between different stellar types, while the slope 
f3 ~ 2.2 in all cases. In case of collisions dominated by 
gravitational focussing, the cross section is of the same 
form, but with slope f3 = 2. Henceforth, we can express 
the tidal capture rate i^id in the same way as the collision 
rate in equation (|6]) , but with closest approach given by 

R min = T(2 - {^]~" Ra, (8) 

7T Ml + M 2 V u *,2 / 

where T is the complete gamma function. 

In what follows, we focus on the collision and tidal cap- 
ture rate between two neutron stars. Since we assume a 
homogeneous distribution of stars in the core, these re- 
sults can be re-scaled for close encounters between other 
types of stellar objects. In Table El we give these scaling 
factors in case of a neutron star (NS), a stellar black hole 
(BH), and a white dwarf (WD), with (typical) masses 
and radii as indicated. After taking into account the rel- 
ative differences in fractional numbers /j, the collision 
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burning. 

Even though it is clear that the det ails of these fluc- 
tuati ons are far from known (but see Heg gie fc Gierszl 
2009), the core after collapse is expected to be on aver- 
age larger than those of the observed post-core-collapse 
GCs, and to increase over time to explain the overall dim- 
mi ng in the central surfa ce brightness (see also Fig. 10 
of IGiersz fc Heggidl2009l ). To mimick this average be- 
haviour, we assume a gradually expanding core after 
deepest core collapse, which in turn we describe by two 
classifcal self-similar solutions in two phases. In the early 
phase expansion mirrors the self-similar late phase co l- 
lapse, but at a slower rate (In agaki fc Lvnden- Bcll 1983). 
After a time since deepest collapse r oughly equal t o the 
duration of the late phase collapse (|Heggid I1985D . the 
core enters the late phase in which it fo llows the self - 
similar expansion of an isolated system ([Henonl fl965). 
We derive the encounter rate during all four phases of 
collapse and expansion. 

In addition to this core evolution model with gradual 
core expansion after deepest collapse, we also consider a 
model in which the core properties are kept fixed after 
deepest collapse. Since the latter "halt of collapse" would 
mimick binary burning without fluctuations, it provides 
an upper limit to the contribution of high-density cores, 
and hence an upper limit to the encounter rate. On the 
other hand, while the core expansion is just an aver- 
age approximation to the possible complex fluctuations 
after deepest collapse, we believe it to provide a closer 
estimate to the encounter rate. Moreover, we find be- 
low (see also Figure [TJ]) that the predicted fraction of 
(post-)core-collapse GCs is similar to that observed for 
the MW. The core expansion could also naturally explain 
why only loose GCs have a shall ow/depleted global mass 
function (|De Marchi et al.ll2007h . 

We adopt M 15 as a proto-typical GC that un- 
derwent core collapse. Because this GC has 
been ext ensively studied and modelled (iDull et all 
19971 iMcNamara. Harrison fc Baumgardti 120041: 
van den Bosch et al.l 120061 ). its current properties 
are known accurately. Given the high concentration 
c = log(r c /r t ) > 3, where r* is the tidal radius, we as- 
sume that the core of M 15 is currently still very close to 
deepest collapse. The fr action of its age i agc ~ 13.2 Gyr 
(jBaumgardt et all I2003D that M 15 spent in the early 
and late phase of core collapse depends on its initial con- 
centration Ci n i. The latter can be inf erred from evo lved 
(single-mass) Fokker-Planck models (jQuinlanl 119961 ) for 
a given initial half- mass relaxation time i rn (0). In turn, 
the latter follows from the current half-mass relaxation 
time i rh ^ 1.42 Gyr as t lh (0) ~ t lh (t agc ) - t ag c/£dis, 
under the assum ption of a constant mas s loss rate 
from the GC fe.g. iVesperini fc Heggid[T997l ). Since the 
rate of dissol ution, £dis, depends again on the initial 
concentration (jGnedin et al.l Il99"9h . we have to find Ci„i 
iteratively. 

Doing this, we infer for M 15 £dis — 51, i r h(0) ~ 
1.68 Gyr, and corresponding c m i ~ 1.70. This implies 
that M 15 spent about 8.0 Gyr in the early phase, and 
took another 5.2 Gyr to reach the current phase of deep 
core collapse. Even though the encounter rate during the 
late phase is expected to dominate over the encounter 
rate during the early phase, we still take the increase in 
v cuc during the early phase into account. In particular, 
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Fig. 13. — An estimate of the encounter rate between two neutron 
stars in M 15 as function of time (or redshift at the top). The 
thin lower curves show the predictions for the collision rates, while 
the thick upper curves present the tidal capture rates. The solid 
curves are for the fiducial core evolution model, whereas the dashed 
curves assumes a faster collapse (and expansion), and the dash- 
dotted curves indicate the effect of an earlier collapse. The dash- 
triple-dotted curve shows an additional model in which the core 
pro perti es are kept fixed after deepest collapse. (See Sections 14.21 
and 14.41 for further details.) 

we take into account the changes in the fractions of stellar 
types due to mass segregation in the early phase, whereas 
we assume the fractions to remain constant during deep 
collapse. We predict the core expansion by mirroring 
the late phase collapse but with an expansion rate that 
is a factor 3 slower than the collapse rate. After about 
5.2 Gyr (the duration of the late phase collapse) in the 
future, we expect the core of M 15 to go into the late 
self-similar expansion phase. 

Figure [13] presents the resulting rates as function of 
time (and redshift along the top axis) for close encoun- 
ters between two neutron stars in the core of M 15. The 
number fraction of neutron stars in the core at the begin- 
ning of the early phase is less than 1%, but it rises due 
to mass segreg ation to about 55 % during the self-similar 
deep collapse ([Dull et all 119971 ). The thin lower curves 
show the predictions for the collision rates, while the up- 
per thick curves present the tidal capture rates. The solid 
curves are for the fiducial core evolution model, whereas 
the dashed curves assumes a faster collapse (and expan- 
sion), and the dash-dotted curves indicate the effect of 
an earlier collapse. These two variations on the fiducial 
model are further discussed in Section FP1 below . Finally, 
the dash-triple-dottcd curve shows the additional model 
in which the core properties are kept fixed after deepest 
collapse. 

4.3. Average encounter rate 

We compute an average encounter rate per host galaxy 
by assuming that all GCs follow an evolution similar 
to that of M 15, with a scaling in time, based on the 

half-mass relaxation time i rn cx rV^M 1 / 2 . We adopt 
the distribution in half- mass radii r^ of GCs derived by 
Uordan et al.l (|2005l ). which is inde pendent from the dis- 
tribution in their total masses M (|McLaughlinll2000D . It 
is believed that initially the GC mass distribution fol- 
lowed a power-law, but that especially the less-massive 
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GCs dis solved, most, if not all, b efore going into core 
collapse ([McLaughlin fc Faiill2008D . Hence, we adopt the 
current GC mass distribu tion, which is well d escribed by 
a lognormal distribution ([Jordan et al.ll2006l ) . After ran- 
domnly drawing from these distributions in r^ and M, we 
arrive at a current half-mass relaxation time, i.e., after a 
time i ago — 13 Gyr since the formation of an average old 
GC. Adopting the same initial concentration c- m i = 1.70 
as estimated for M 15, we infer the corresponding initial 
half-mass relaxation time, which we use to re-scale the 
above derived encounter rate evolution for M 15. 

The current and initial half-mass relaxation distribu- 
tions are shown as solid and dotted histograms in the 
left panel of Figure Qj] while the vertical lines indicate 
the corresponding values for M 15. Since M 15 has a 
larger than average half-mass relaxation time, we expect 
that currently more than half of all GCs are already past 
reaching deepest core collapse. This is confirmed by the 
thin solid curve in the right panel of Figure 1141 showing 
the fraction of GCs that past reaching deepest core col- 
lapse as function of time. In this fiducial core evolution 
model it is likely that some time after deepest core col- 
lapse the GCs are not anymore observed to be post-core- 
collapse because the core is expanding. Henceforth, to 
estimate the fraction of observed post-core-collapse GCs, 
we only count at a given time GCs that are in the late col- 
lapse or early expansion phase when the core is smaller. 
The resulting thick curves show that for the first three 
models the current fraction (vertical solid line) is simi- 
lar to the fraction of ~ 20 % of post-core-collapse GCs 
observed in the Milky Way. Only the halt-of-collapse 
model predicts a much higher fraction because the core 
properties are held fixed after deepest collapse. 

The number of GCs for which we re-scale the encounter 
rate evolution of M 15 depends st rongly on the tota l lu- 
minosity Lgai of the host galaxy ([McLaug hlin 199^). In 
Figure [TS] we plot the combined encounter rate for a typi- 
cal L-k galaxy with about a thousand GCs, as well as for a 
galaxy that is a factor 5 less and more luminous with the 
number of GCs comparable to that in the Milky Way and 
in a giant elliptical galaxy, respectively. In the latter two 
cases we show only the tidal capture rates for the fidu- 
cial model, wheras for the L-k galaxy we also indicate the 
effect of a faster and earlier collapse as discussed in Sec- 
tion 14.41 as well as the additional half-of-collapse model. 
The results for the collision rate are not shown for clar- 
ity, but are similar except for a factor ~ 6.7 decrease 
in amplitude. In all cases, the encounter rate increases 
strongly up to redshift z ~ 1, but slowly decreases again 
toward z = 0. 

Next, we estimate a global encounter rate by com- 
bining the previously derived encounter rates per host 
galaxy with the galaxy luminosity density distribution. 
The latter is well described by the Schechter function 

$ ga i(L) = ($ gal * I L-k) (L/L*) a exp(-L/L*), (9) 

with $ ya1 * = 3.77 x 10~ 3 Mpc- 3 and a = -1.30 
([Faber et al.ll200"l . Since the above estimate of the com- 
bined encounter rate evolution depends on the current 
luminosity of the host galaxy, the evolution of & ga i(L) 
is not needed to compute the global encounter rate as 
function of time. However, the latter evolution should 
be taken into account when discussing the properties of 



host galaxies, which do change significantly with redshift 
(see Section [Q)l . 

In Figure 1161 we present our prediction for the global 
encounter rate -R en c (in yr -1 Gpc -3 ) between two neu- 
tron stars in the core of GCs. The rates both of collisions 
(thin curves) and of tidal captures (thick curves) show 
a clear increase around z ~ 1, followed by a gradual 
decline. The tidal capture rate for the fiducial model 
(thick solid curve) peaks around z ~ 0.7, at a value 
of i?enc — 55 yr ~ 1 Gpc -3 . This value is fully consis- 
tent with the estimated event rates for SRGBs, which 
are of the order 8-30 y r" 1 Gpc~ 3 for isotropic radiation 
([Guetta fe Piran 20Q6|) . Collimated radiation with open- 
ing angles of 30-60°, would boost up the event rates of 
SGRBs to our predicted global encounter rates. How- 
ever, since the latter are likely to be rather conservative, 
smaller opening angles would still be possible. 

4.4. Faster and ealier collapse 

To infer the initial concentration of Ml 5 and the cor- 
responding time spent in the early and late phase of col- 
lapse, we used the res ults of a series of cluster models 
with single- mass stars (Quinlan I1996I ). Several compli- 
cating factors in real clust ers can change the rate of col- 
lapse (see e.g. review by iChernofil [1993): for example 
stellar evolution and especially the presence of primor- 
dial binaries slow down the collapse, while mass loss from 
a tidal boundary and in particular mass segregation ac- 
celerate the collapse. Even though the fraction of pri- 
mordial binaries in GCs is believed to be only a few per 
cent fe.g. lHurlev et al.ll2007Hbavis et al.ll2008l) . the bind- 
ing energy of a single hard binary is sufficient to signifi- 
cantly slow down or perhaps even prevent core collapse 
([Hut et al.lll99llHeggie et al.ll2006f ). On the other hand, 
because the relaxation time is inversely proportional to 
the mean stellar mass, the acceleration of core collapse 
due to mass segragation is roughly proportional to the 
increase in the mean mass, w hich in the core of M 15 is 
about a factor 1.19/0.43 ~ 3 ([Dull et al.lll99l . 

The detailed inclusion of these competing mechanism 
is beyond the scope of this paper, but to investigate the 
effect of a changing collapse rate on the encounter rate, 
we assume a 50% faster core collapse (and expansion) 
than the fiducial M 15 core evolution model. At the 
same time, this implies a lower initial concentration of 
Ci n i ~ 1.47. The results of this "faster collapse" model 
on the encounter rate are shown with dashed curves in 
Figures HBTHBj and [TH As can be seen from Figure [TH 
the main effect is a lowering of the global encounter 
rate by about 24% with respect to the fiducial model 
to -Rene — 41yr~ 1 Gpc -3 at the peak, which remains 
around z ~ 0.7. 

We already mentioned in Section 14.21 that after deep- 
est core collapse the core can un dergo fluctuations, either 
due to gravothermal o scillations ([Bettwieser fc Sugimotol 
[IMilGoodma^lT98l . or as a result of the stochastic na- 
ture of binary burning ([Giersz fc Heg gic 2009). In either 
case, these fluctuations can result in change in the core 
radius r c by factors of a few, and even order of magni- 
tude changes in the core density p c . We do not attempt 
to capture this fluctuating nature, but we consider here 
a possible consequence of it. Because the (central) re- 
laxation time is inversely proportional to p c , only a rel- 
atively small fraction of the time during an oscillation 
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Fig. 14. — The left panel shows the distribution of half-mass relaxation times t T h <x r^J 2 M 1 / 2 , based on the independent distributions 
in half-mass radii and total mass M of globular clusters. This results in the solid histogram, which turns into the dotted histogram 
for initial half-mass relaxation times t r h(0). The vertical lines show that M 15 has a larger than average half-mass relaxation time, so we 
expect that currently more than half of all GCs are already past reaching deepest core collapse. This is confirmed by the thin solid curve 
in the right panel, showing the fraction of GCs that past reaching deepest core collapse as function of time. However, in this fiducial core 
evolution model it is likely that some time after deepest core collapse the GCs are not anymore observed to be post-core-collapse because 
the core is expanding. Henceforth, to estimate the fraction of observed post-core-collapse GCs, we only count at a given time GCs that 
are in the late collapse or early expansion phase when the core is smaller. The resulting thick curves show that for the first three models 
the current fraction (vertical solid line) is similar to the fraction of ~ 20 % of post-core-collapse GCs observed in the Milky Way. Only the 
fourth halt-of-collapse model predicts a much higher fraction because the core properties are held fixed after deepest collapse. 
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Fig. 15. — An average encounter rate between two neutron stars 
in the core of globular clusters per host galaxy with luminosity 
L ga l = {0.2, 1.0, 5.0} X L-k. The expected number of GCs per host 
galaxy, given in brackets, follows from McLaughlin (1999). For 
each globular cluster, the collision rate of M 15 is re-scaled by the 

half-mass relaxation time t r ^ oc r^ 2 M 1 / 2 , randomly drawn from 
the independent distributions in half-mass radii and total mass 
M of globular clu sters. The meaning of different line styles is the 
same as in Fig|13l 



Fig. 16. — The predicted global rate of close encounters between 
two neutron stars in the cores of globular clusters, per volume of 
space and as function of time (or redshift at the top). The values 
are compute d fr om combining the average collision rate per host 
galaxy (Fig. 1156 with the galaxy luminosity density dis trib ution. 
The meaning of different line styles is the same as in Fig |13l 
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will be spent at large p c and correspondingly small r c . 
Still, if we happen to observe a post-collapse GC dur- 
ing this special time at the "peak" of the fluctuation, 
we will underestimate the time elapsed since deepest col- 
lapse under our assumption of a steady expansion in our 
fiducial core evolution model. To investigate the possi- 
ble effect on the encounter rate, we consider a shift the 
time of deep core collapse back in time by an amount 
of 2.6 Gyr, i.e., half of the duration of the late-phase 
collapse. As a consequence, the early-phase collapse is 
shorter by the same amount and the initial concentration 
increases to c m i ~ 1.79. The smooth expansion model in 
this case predicts a current core radius and core den- 
sity that are respectively almost 65% larger and a fac- 
tor 3 smaller than the ob served values of r c = 0.05 pc 
(iNovola fc Gebhardtl|2006D and p c = 6.5 x 10 6 M© pc" 3 
(jvan den Bosch et al.ll200a ). The results of this "earlier 
collapse" model on the encounter rate are indicated by 
the dash-dotted curves in Figures ITBTTlBI and H21 Fig- 
ure [16] shows that the global encounter rate is increased 
by about 14% with respect to the fiducial model to 
i? cnc ~ 62yr _1 Gpc~ 3 at the peak, which occurs ear- 
lier, at z ~ 1.0. The latter is directly related to the 
applied shift in time, whereas the former is the result of 
more GCs reaching deep collapse in i age = 13 Gyr (see 
also Figure [T4"|) . 

The above faster and earlier collapse models imply re- 
spectively a lower and higher initial concentration for 
M 15 than used in the fiducial model. For each model, 
we then fixed the initial concentration of all GCs to that 
of M 15 when computing the average encounter rate. It is 
already very difficult to infer the initial concentration for 
M 15, let alone the distribution of initial concentrations 
for the population of old GCs. Even so, we have recom- 
puted the average encounter rate adopting a Gaussian 
distribution in c m i 1 with a mean at the initial concentra- 
tion of M 15 and with dispersion in the range [0, 0.3]. As 
a result of the dispersion, the rate of dissolution £dis is 
not anymore fixed to that of M 15, which in turn leads 
to an additional broadening in the half-mass relaxation 
distribution when converting from current t r h to initial 
irh(0) values. However, since the resulting variations in 
the conversion are typically less than ~ 0.1 Gyr, the ef- 
fect on the average encounter rate is very small. Chang- 
ing the mean, as we do in the faster and earlier collapse 
models, does have a significant effect on the average en- 
counter rate. Nevertheless, we see from Figure [16] that 
the evolution of the global encounter rate is quite similar 
even for significant deviations from the fiducial mdoel. 
In other words, as long as M 15 is indeed representative 
of an (old) GC that underwent core collapse and expan- 
sion, we expect our predictions for the encounter rates 
to be robust. 

4.5. Binaries 

We have not included the potentially important effect 
of binaries in our encounter rate calculations. The frac- 
tion of primor dial binaries in GC s is thought to be only 
a few per cent (jHurlev et al.ll2007j ) , as confirmed by mea- 
suring p hotometric binari es beyond the half-light radius 
(see e.g. lDavis et al.|[2~008[ ). In the core of GCs tidal cap- 
tures can significantly increase the fraction of binaries. 
Placing useful observational constraints on the fraction 
of binaries in the core is not an easy task because many 



of them could be dark as a result of binary-single inter- 
actions in which the more numerous main sequence stars 
(MS) and white dwarfs (WD) are replaced by compact 
remnants. Binary systems have a larger cross section 
for intersection and they mass segregate faster into the 
denser core than individual stellar objects. Three-body 
interactions are thus indeed likely to create compact bi- 
naries in the cores of GCs, including double neutron star 
(NS-NS) and black hole-neutron star (BH-NS) systems. 
As already noted in the Introduction, the merger of such 
dynamically formed compact binaries is potentially an 
important ch annel for the producti on of SGRB progen- 
itor systems (jGrindlav et al.l I2006D , different from the 
close encounters between single compact stellar o bjects 
we have considered so far. As we discuss below in 14.5.21 
for encounters involving the exchange of a main sequence 
star with a neutron star or a black hole, the resulting sep- 
aration will be too large for coalescence to take place in 
less than a Hubble time. Thus this channel is unlikely to 
contribute to the overall formation rate of viable SGRB 
progenitors. The occurrence of a direct impact may re- 
lease a substantial amount of energy, but whether this 
can lead to the conditions required for SGRB produc- 
tion remains to be evaluated. 

The three main ways to arrive at compact binaries 
are (i) primordial compact binaries, (ii) primoridal bi- 
naries of which the non-compact member(s) are replaced 
by compact objects via one (or two) exchange interac- 
tions, and (iii) the same as (ii) but starting from binaries 
which themselves are first formed through tidal capture 
of two single stellar objects. We consider all three ways 
below and show that in the cores of GCs, the latter dom- 
inates the formation rate of compact binaries, but is not, 
as commonly believed, dominant over the rate of close 
encounters between two single compact objects. While 
SGRBs are potentially created instanteneously when two 
compact objects closely encounter each other, the sepa- 
ration of the compact binary has to be small enough such 
that the time for the compact objects to merge, in addi- 
tion to the formation time, is less than the Hubble time. 

Note that we expect very little contribution to the close 
encounter rate from interactions between a single com- 
pact object (NS or BH) and a compact object that is a 
member of a binary which also contains a extended ob- 
ject (MS or WD). Even though we show below that the 
fraction of the latter binaries can be(come) significant, in 
nearly all cases the (exc hange) interaction will be wi th 
the extended object (e.g. ISigurdsson fc Phinnevlll993l ). 

4.5.1. Primordial compact binaries 

The formation of primordial binaries follows the cosmic 
star formation rate which pea ks at high redshift z ~ 3 
(e.g. lHopkins fc Beacomir2006| ) . Primordial binaries with 
two massive stars (> 8M Q ) that survive the mass loss 
and possible asymmetry of the supernova explosion of the 
secondary provide compact NS-NS and BH-NS binaries. 
They potentially could produce SGRBs at high redshift 
with a distribution extending to lower redshift depending 
on the range of initial binary separations and thus merger 
times. However, even if a non-neglible fraction of primor- 
dial massive binarie s survive and result i n close compact 
binaries (~ 10~ 2 5 . iNaravan et al.l [199 11 . massive bina- 
ries are only a very small fraction of all primordial bi- 
naries for a reasonable, not too top-heavy initial mass 
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function (~ 10 -4 5 , for an IMF with a Salpeter power 
law slope from 100 M Q down to 0.3 M©, and constant 
down to 0.1 M©). In GCs with already at most a few per 
cent primordial b inaries (~ 10~ 1,5 . IHurlev et ail 120071 : 
iDavis et al.ll2008l) and likely additional co mplications to 
form and retain compact binaries (see also llvanova et ahl 
2008), their possible contribution to producing SGRBs 
should be negligible; the GCs in the Milky Way contain 
a total of ~ 10 7 ' 5 stars, which implies that we expect 
only ~ 1 primordial compact binary formed in the GCs 
of the Milky Way. 

However, in the "field" of a galaxy the primordial bi- 
nary fraction is of order unity and the vast amount of 
them results in a number of compact binaries that is 
not anymore insignificant; we expect ~ 10 5 primordial 
compact binaries formed in th e "field" of the Milky Way 
(see, e.g. iNaravan et ahl 119911 : lPhinnevlll99lL for initial 
estimates). Merger rate calculations for such primordial 
binaries in t he galactic field have become increasingly so- 
phist icated ( Kalog era et al.1 l2004t IQ'Shaughnessy et al.1 
120051 [20081 12010D . and now consider not only various stel- 
lar evolutionary channels, but also different host galaxy 
types (spiral vs. elliptical). The results are still hin- 
dered by the fact that there is a small number of ob- 
served systems, giving broad estimates in the range 1— 
lOOyr^Gpc" 3 . 

4.5.2. Dynamically formed compact binaries 

Compact binaries formed dynamically in GCs have 
long delay times between compact object birth and bi- 
nary formation as they ha ve to wait for collapse of the 
core (|Hopman et al.l [20061) . This helps their formation 
in three ways: (i) mass segregation increases the relative 
fraction of heavier-than-average binaries as well as com- 
pact remnants in the core, (ii) tidal capture increases the 
binary fraction in the core, in particular the "hard" bina- 
ries with small separations (whereas the "soft" binaries 
with large separations are being "ionized" ) , (iii) the high 
core density allows three-body exchange interactions in 
which a member of an existing binary is replaced with 
a compact remnant. Henceforth, we expect in the cores 
of GCs the compact binary formation rate to follow the 
close encounter rate, with a similar peak around z ~ 0.7 
as shown in Figure [16J Whereas after this initial delay 
the potential production of SGRBs is nearly instanta- 
neous in the proposed encounter scenario, there is an ad- 
ditional delay in the binary merger scenario given by the 
rate of angular momentum losses through gravitational 
waves. 

For a NS-NS binary dynamically formed around z ~ 
0.7 to merge within the next ~ 6 Gyr through the emis- 
sion of gravitational waves, the binary separation re- 
quired is a < 3.93 R©. As mentioned above, most of 
the compact binaries in GCs are expected to be formed 
through three-body exchange interactions: a NS (average 
mass 1.4 M©, and average radius 10 km) replaces the less- 
massive main-sequence star (MS; 0.4M©, 3.3 x 10 5 km) 
or white dwarf (WD; 0.5 M©, 1.1 x 10 4 km) member in 
an existing binary with a NS. The latter NS-MS or NS- 
WD binary might be itself the result of a previous simi- 
lar three-body exchange interaction with a MS-MS, WD- 
MS or WD- WD binary. We assume as minimum separa- 
tion of these binaries the closest distance before ignition 



of mass transfer occurs (|Paczvhskiiri97lD . This yields 
a > 1.72R© and a > 0.054R© for NS-MS and NS-WD 
binaries, respectively. We obtain similar lower limits of 
about 1.312, 1.365, and 0.043 in units of R© for MS-MS, 
WD-MS, and WD- WD binaries, respectively. 

Even though after each exchange the binary hardens, 
the binary separation after an exchange, ag n , is typically 
large r than the binary separation before the exchange, 
aini- IBigurdsson & Phinncy (1993) show that the distri- 
bution of Ofln/tiini strongly peaks around the mass ratio 
mf/m,e, between the single stellar object from the "field" 
and the stellar object "exchanged" from the binary, but 
with a long tail towards lower values. Taking for aj n i the 
minimum separation of each binary, we compute the me- 
dian of the latter distribution^ to arrive at the average 
ag n after the first exchange interaction, and repeat this 
for the MS-MS, WD-MS, and WD- WD binaries that re- 
quire a second exchange interaction to arrive at a NS-NS 
binary. For all three binaries with a MS as a member the 
resulting minimum binary separations (5.21, 12.8 and 
12.4 in units of R© for NS-MS, MS-MS, and WD-MS, 
respectively) are larger than the above maximum binary 
separation for merger of the dynamically formed NS-NS 
binary (3.93 R© for merging within ~ 6 Gyr). Note that 
this is still the case even when we ignore the time to 
reach core collapse and allow the merger to take a full 
Hubble time, corresponding to a maximum binary sepa- 
ration of 4.8 R©. On the other hand, a NS-WD binary 
results in a NS-NS binary with a minimum separation 
a > 0.126 R© well below the maximum for merging. Sim- 
ilarly, a WD- WD binary yields a NS-WD binary with 
a > 0.102 R© after the first exchange, and a NS-NS bi- 
nary with a > 0.240 R© after the second exchange. 

We assume that the formation rate of these binaries 
through three-body exchange interactions follows that of 
the collision rate, but with a different amplitude, which 
we estimate below. 

4.5.3. Compact binaries from primordial binaries 

Here we estimate the rates at which compact bina- 
ries in the cores of GCs are dynamically formed from 
primordial binaries, of which one or both members are 
non-compact stellar objects. We take / pr i = 0.02 for 
the total fraction of primordial binaries in GCs, which 
is not wel l known but thought to be at most a few per 
cent (e.g. IHurlev etal] [20071 [Davis et all [2001 . To de- 
rive the fraction of primordial binaries that have a NS 
and a WD, we adopt an IMF with Salpeter power-law 
slope from 100 M© down to 0.3 M©, and constant down 
to 0.1 M©. Supposing that stars above 20 M© evolve into 
black holes (BHs), stars above 7M© turn into SNe, stars 
above 0.8 M© end their life as WDs, and lower-mass stars 
are still in the MS, we find corresponding number frac- 
tions of about 0.0012 0.0043 0.1014, and 0.8930 for BH, 
NS, WD, and MS stars, respectively. While we start 
with drawing the masses of the primary and secondary, 
Mi and M2, indepdently from the IMF, we also con- 
sider the case that the mass ratio M2/M1 is distributed 
uniformly in the range from zero to unity. The former 
case is fully consistent with observed (spectroscopic) bi- 

3 We adopt an approaching velocity from infinity of = 2tr c ~ 
20 km s , but the results are not sensitive to the precise value that 
is adopted. 
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naries with longer periods (P > 10 3 days), but there 
are indications of a flatter distribution in mass ratios for 
the shorter-period binaries, possibly because stars that 
form close to each other become closer in mass by reach- 
ing an equilibrium in mass transfer. We adopt a log- 
normal distribution in binary periods P (in days) with 
mean ii\ os p = 4.8 and stan dard deviation ai og p = 2.3 
puquennov fc Mavorj["l991h . 

We first consider the formation of NS-NS binaries from 
primoridal NS-WD binaries through a single exchange 
interaction. With the masses of the primary and sec- 
ondary independently drawn from the IMF, we expect 
a fraction /NS-WD.pri — 8.81 x 10~ 4 of all primordial 
binaries to be a NS-WD binary. Together with the 
stars and compact remnants that are heavier than the 
MS stars, the binary sinks to the core of a GC. Dur- 
ing deep core collapse the NS-WD undergoes a three- 
body exchange interaction with a NS for whic h we adopt 
the n umber fraction / NS ~ 0.547 of M 15 pull et al.l 
1997). We only count the NS-WD binaries with sepa- 
rations in the range [0.054, 1.669] R©. The lower limit 
is set to avoid stable mass transfer, while the upper 
limit increases after the exchange interaction to the 
maximum separation of 3.93 R© allowed for merger of 
the dynamically formed NS-NS binary within ~ 6 Gyr. 
Given the above lognormal distribution in binary peri- 
ods, the fraction of NS-WD binaries in this range of 
separations is / scp — 1.83 x 10~ 2 . The corresponding 
mean separation is a ~ 0.829 R©. We now obtain the 
rate from eq. © with i? m j n = a + tns — 0.829 R©, 

/l = /pri /NS-WD,pri /scp — 3.23 X 10~ 7 , f% = /nS — 

0.547, and mi + m 2 ~ 2 x 1.4 + 0.5 = 3.3 M©. In 
this way, we predict that the rate i?dyn to turn a pri- 
mordial NS-WD binary dynamically into a NS-NS bi- 
nary that can merge within the next ~ 6 Gyr, is a factor 
Rdyn/Rco\ — 4.02 x 10~ 2 smaller than the collision rate, 
and a factor i?dyn/-Rtid — 0.618 x 10~ 2 smaller than the 
tidal capture rate. 

We now turn to the additional dynamical formation 
of a NS-NS binary from a primoridial WD- WD binary 
through two subsequent three-body exchange interac- 
tions with a NS. We assume that the first exchange 
to form the NS-WD binary happens some time before 
deepest core collapse, with a formation rate that can 
be computed in the same way as above. The fraction 
of primodial binaries that are WD- WD binaries is about 
/wD-WD.pri ^ 1-03 x 10~ 2 . We count the WD- WD bina- 
ries with separations in the range [0.043, 0.708] R©, with 
the lower limit again to avoid stable mass transfer, while 
the upper limit increases after two subsequent exchange 
interactions to the maximum separation of 3.93 R© al- 
lowed for merger of the dynamically formed NS-NS bi- 
nary within ~ 6 Gyr. This range in separation corre- 
sponds to a fraction / scp ~ 1.10 x 10 -2 and mean sepa- 
ration a ~ 0.371 R©. We then use as before eq. (HJ) with 
R miD ~ 0.371 R Q , while we set f x = / pri /ns-wd / scp - 
2.27 x 10~ 6 , and m 1 +m 2 -1.4+2x 0.5 = 2.4M©. As 
in Section 14.21 we assume f% — /ns (t) to vary from its 
primordial value to 0.547 during the early-phase collapse, 
and then to remain the same during the late-phase col- 
lapse. We intergrate the resulting rate of exchange inter- 
actions between a primordial WD- WD binary and single 
NS over the full duration of the core collapse phase. This 



yields the expected number (~ 4.24 x 10~ 4 ) of dynam- 
ically formed NS-WD binaries with separations within 
the above range. We then divide this number by the 
number of stellar systems (~ 2.86 x 10 3 ) in the core 
N c = (47r/3)n c r 3 at deepest core collapse to arrive at 
a fraction /NS-WD,dyn — 1-48 x 10~ 7 . Next, we repeat 
the rate computation but now for dynamically forming 
a NS-NS binary in a second exchange interaction with 
h = /ns-wd, dyn, ji = /ns - 0.547, ra x + m 2 ^ 2.4 M©, 
and i? m in = 0.873 R©. The latter value arises from the 
above mean separation a ~ 0.371 R© after the first ex- 
change interaction, multiplied by the (median) increase 
after the second exchange interaction. In this way, we 
find that the rate to dynamically form a NS-NS binary 
from a primordial WD- WD binary through two subse- 
quent three-body exchange interactions is smaller than 
the collision rate by a factor i?dyn/-Rcoi — 1-94 x 10~ 2 , 
and smaller than the tidal capture rate by a factor 

i?dyn/-Rtid ^ 0.299 X 10- 2 . 

By combining the above two formation channels, we 
thus expect that the rate to dynamically turn primordial 
binaries into compact NS-NS binaries that can merge 
in time is smaller than the tidal capture rate by a factor 
~ 0.916 x 10~ 2 . In case that the mass ratio M 2 /Mi of the 
primordial binary is uniformly distributed, the fraction 
of primordial NS-WD and WD- WD binaries increases to 
/NS-WD,pri - 4.76 xl0~ 3 and / W D-WD, P ri ^ 7.18 xl0~ 2 , 
i.e., a factor ~ 5.41 and ~ 6.98 higher than when Mi 
and M2 are independently drawn. As a result, the rate 
to form compact NS-NS binaries increases by a factor 
~ 5.92, but is still only a fraction ~ 5.43 x 10~ 2 of the 
tidal capture rate. 

A similar calculation for compact BH-NS binaries dy- 
namically formed from primordial binaries shows that 
also in this case the minimum separations of primordial 
binaries with a MS star as a member are larger than the 
maximum separation of < 6.34 R© for which the BH-NS 
can merge within ~ 6 Gyr. However, BH-NS binaries can 
form dynamically from primordial NS-WD or BH-WD 
binaries via a single exchange interaction with respec- 
tively a BH or NS, or from primoridal WD- WD binaries 
via two subsequent exchanges interaction with both a 
NS and a BH. We compute the rate for each of these 
formation channels in the same way as above. We as- 
sume that the fraction of BHs relative to the (changing) 
fraction of NSs remains equal to the primoridal fraction, 
i-e., /bh//ns — 0.283 throughout. If instead of chos- 
ing the masses of the primary Mi and secondary M 2 of 
the primordial binary indepedently, we assume that the 
mass ratio M2/M1 is uniformly distributed, the fraction 
of primordial NS-WD, BH-WD and WD- WD binaries in- 
creases by a factor 5.41, 1.73, and 6.98, respectively. At 
the end, the combination of the three formation chan- 
nels yields a formation rate of compact BH-NS binaries 
is smaller than the tidal capture rate between a single 
BH and NS by a factor ~ 0.593 x 10~ 2 if Mi and M 2 are 
independently drawn, and a factor ~ 1.88 x 10~ 2 when 
Mil Mi is uniformly distributed. 

4.5.4. Compact binaries from tidally captured binaries 

It is clear from the above calculations that in the cores 
of GCs close encounters between single compact objects 
are significantly more likely than mergers of compact bi- 
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naries that are dynamically formed from primordial bina- 
ries. When starting instead from binaries that are formed 
through tidal captures the chance of dynamically form- 
ing and merging a compact binary can be much higher 
as we show below. As before binaries with a MS as a 
member are omitted, because they result after exchange 
interaction(s) in compact binaries with minimum sep- 
arations that are larger than the maximum separation 
allowed for them to merge in time. 

We first consider the formation of NS-NS binaries from 
tidally captured NS-WD binaries. To compute the frac- 
tion of tidally captured NS-WD binaries, we integrate 
the tidal capture rate between a NS and WD over time 
until deepest core collapse, and then divide by the num- 
ber of stellar systems in the core. The tidal capture rate 
is given by eq. (j6|) after substitution of R m i n from eq. (|8]l . 
We assume fi = /ns(£) and fi — fwu(t) to vary dur- 
ing the early-phase collapse from their primordial values 
to respectively 0.547 and 0.218, and then to remain the 
same during the late-phase collapse. This yields the ex- 
pected number (~ 8.86) of tidally captured NS-WD bi- 
naries with (pericenter) seperations up to O.1O9R0. This 
is well below the upper limit of 1.669 R© that increases 
after the exchange interaction to the maximum separa- 
tion of 3.93 R© allowed for merger of the dynamically 
formed NS-NS binary within ~ 6 Gyr. The lower limit 
to avoid stable mass transfer is 0.054 R©, and gives rise 
to a fraction / scp ~ 0.511 of tidally captured NS-WD 
binaries with separations in this allowed range, and with 
corresponding mean separation of a ~ 0.093 R©. Finally, 
we divide by the ~ 2.86 x 10 3 stellar systems in the core 
to arrive at a fraction /NS-WD,tid — 1-72 x 10~ 3 of tidally 
captured NS-WD binaries. The calculation for the subse- 
quent exchange reaction with a NS, is the same as above 
in Section l4.5 3l for the formation of NS-NS binaries from 
primordial NS-WD binaries, except that in eq. (|6|) we 
now have /1 = /NS-WD.tid- In this way, we estimate 
that the rate i?dyn to turn a tidally captured NS-WD 
binary dynamically into a NS-NS binary that can merge 
within the next ~ 6 Gyr is a factor i?dyn/-Rtid — 3.70 
larger than the tidal capture rate. 

We now turn to the additional dynamical formation 
of a NS-NS binary from a tidally captured WD- WD bi- 
nary through two subsequent three-body exchange in- 
teractions with a NS. Since the tidal capture rates be- 
tween two WDs and between a WD and a NS are simi- 
lar (Table [3]), we expect the rate to be significantly less 
than just estimated for NS-NS binary formation from a 
tidally captured NS-WD binary, which only needs one 
subsequent exchange interaction. We assume that both 
the tidal capture of the two WDs and the first sub- 
sequent exchange happen before deepest core collapse. 
This means we first calculate the fraction of WD- WD bi- 
naries formed through tidal capture as function of time, 
i.e., /wD-WD,tid(i)- The calculation is the same as for 
/NS-WD,tid above, but with fi = f 2 = fwo(t), and we 
integrate the tidal capture rate up to time t, followed by 
division with the number of stellar systems in the core 
N c (t) = (4Tr/3)n c (t)r c (t) 3 at that time. To have the 
WD- WD binary (pericenter) seperations up to 0.092 R© 
within the allowed range of [0.043, 0.708] R©, implies 
/ sop ~ 0.530 and a mean separation of a ~ 0.072 R©. 
Next, we substitute i? m in = a + tns — 0.072 R©, /1 = 
/wD-WD,tid(0; /2 = hs{t) andmi+m 2 ~ 2x0.5+1.4 = 



2.4 M© in eq. ([5]). Integrating the resulting rate over the 
full duration of the core collapse phase, yields the ex- 
pected number (~ 4.08 x 10 -3 ) of dynamically formed 
NS-WD binaries from tidally captured WD- WD bina- 
ries. We then divide by the ~ 2.86 x 10 3 stellar systems 
in the core at deepest collapse, to arrive at a fraction 
/NS-WD,dyn — 1-43 x 10~ 6 . As above in Section 14.5.31 
we repeat the rate calculation for dynamically forming 
a NS-NS binary in a second exchange interaction with a 
NS. In this way, we find that the rate to form a NS-NS 
binary dynamically from a tidally captured WD- WD bi- 
nary through two subsequent three-body exchange inter- 
actions is smaller than the tidal capture rate by a factor 
RdyJRm ^ 0.555 x 10~ 2 . 

By combining the above two formation channels, we 
thus expect that the rate to dynamically turn tidally 
captured binaries into compact NS-NS binaries that can 
merge in time is larger than the tidal capture rate by 
a factor ~ 3.71. As expected, this factor is mainly the 
result of the first formation channel, of which the rate is 
much larger, by a factor ~ 667, than the second forma- 
tion channel. Next, we estimate in a similar way also the 
rates to dynamically form BH-NS binaries from tidally 
captured NS-WD or BH-WD binaries via a single ex- 
change interaction with respectively a BH or NS, or from 
tidally captured WD- WD binaries via two subsequent ex- 
changes interaction with both a NS and a BH. From the 
three formation channels together, we predict a forma- 
tion rate of BH-NS binaries that is larger than the tidal 
capture rate between a single BH and NS by a factor 
~ 8.42. Again, the contribution from the third formation 
channel that involves two subsequent exchange interac- 
tions is much smaller, by factors ~ 72.8 and ~ 371, than 
the first and second formation channel. 

4.5.5. Merger versus encounter rate 

The above calculations provide an estimate of the rela- 
tive frequencies of (binary) mergers and close encounters 
between two compact objects in the cores of GCs. Com- 
bining the results from Sections l4.5.1hl4.5.41 we estimate 
for a compact NS-NS (BH-NS) binary a higher merger 
than tidal capture rate by a factor ~ 3.7 (~ 8.4). 

These estimates can be further improved by allowing 
a range in WD and MS star masses (instead of adopt- 
ing mean masses of 0.5 M© and 0.4M©, respectively), by 
sampling the full distribution in the ratio of binary sep- 
arations before and after a three-body exchange inter- 
action (instead of taking the median value of afi n /ai n i), 
and by including the formation and subsequent merger 
of compact binaries throughout the lifetime of a GC in- 
stead of concentrating on deepest core collapse with a 
peak around z ~ 0.7). Also the binaries will interact 
among themselves, so that four-body interactions could 
become important, especially since the number of tidally 
captured binaries is increasing; when for all ten possible 
binary pairings from a BH, NS, WD and MS, we inte- 
grate the corresponding tidal capture rate up to the time 
of deepest collapse, we find a total of ~ 207 dynamically 
formed binaries. Given the ~ 2.86 x 10 3 stellar systems 
in the core, we obtain a binary fraction of ~ 7.23 % com- 
pared to a primordial binary fraction of a few per cent. 
These and other improvements are clearly interesting but 
beyond the scope of this paper. 

From the above estimates, we can conclude that in the 
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cores of GCs the rate of close encounters between two 
single compact objects can become similar to the rate of 
merging of two compact objects in a binary. The lat- 
ter compact binaries most likely are formed through a 
three-body exchange interaction with a binary which it- 
self originated from a tidal capture between a compact 
and non-compact (WD) stellar object (Section I4.5.4[) . 
The alternative dynamical formation of compact binaries 
from primordial binaries (Section 14.5. 3|) is less likely by 
about two orders of magnitudes. A main reason is that 
the fraction of primordial binaries in GCs is at most a 
few per cent, which in turn also implies that the number 
of primordial compact binaries in GCs is negligble, while 
they are thought to be the dominant source of SGRBs in 
the "field" of a galaxy (Section EH}. while SGRBs are 
commonly assumed to orginated from compact binary 
mergers, we thus find that close encounters between two 
compact objects in the core of GCs can also provide a 
significant contribution. 

5. DISCUSSION 

Dense stellar systems in some ways are like an ecologi- 
cal network, where feedback is extremely important and 
apparently isolated events can have far-reaching conse- 
quences. In a globular cluster, as we discussed, dynami- 
cal interactions between passing stars can form new bina- 
ries and modify the properties and even the membership 
of existing binaries. Motivated by this, we have inves- 
tigated the production of compact binaries via two and 
three body encounters. We find that event rates within 
globular clusters are expected to be significant, and can 
become similar to the overall production of merging com- 
pact binaries. This hints at the underlying possibility 
that SGRB progenitors may not be entirely restricted to 
the most widely favored scenario involving the merger of 
compact binaries in the field. Much of our effort in this 
section will therefore be dedicated to determining what 
are the expected characteristics of SGRBs arising from 
these encounters and how do they compare to recent ob- 
servational constraints. 

5.1. Prospects for the production of SGRBs 

It is clear from the calculations presented in Section l2~4l 
that the formation of an accretion disk around the pri- 
mary is a robust result, regardless of the initial orbital 
parameters (as long as disruption occurs, of course). The 
particulars of each case, however, are variable in several 
aspects, and can lead to interesting diversity. 

First, and most importantly, the resulting disk mass is 
not always equal. This is crucial in terms of the avail- 
able energy in the system, as it sets the overall scale for 
an accretion powered event that could produce a SGRB. 
We find that even for the relatively small variation in 
angular momentum (or equivalently, impact parameter) 
for black hole-neutron star encounters the resulting disk 
mass, Mdi s k, varies almo st by a factor of three. We 
have previously estimated (ILee. R amirez-Ruiz fc Granotl 
120051 : ILee fc Ramirez-Ruizi 120071) that the energy that 
can be potentially extracted from such disks to produce 
a GRB scales as -M| isk for neutrino powered events, and 
as -Mdisk for magnetically dominated bursts. The range 
in mass thus possibly spans an order of magnitude when 
converted to total energy release. 



Second, the nature of the primary itself can produce 
a different outcome. The "cleanest" scenario involving a 
black hole leads to accretion in the hypercritical regime 
discussed above. However, for double neutron star colli- 
sions, a range of possibilities remains. For the one cal- 
culation of this type we have performed, at the end of 
the simulation the central core is rapidly rotating and 
is surrounded by a massive envelope. If the core can 
avoid collapse, it is possible that the rapid rotation will 
wind up the magnetic fi eld to large magnetar-like values 
([Price fc R osswog 2006) a nd all o w for repeated episodes 
of en ergy release ( jUsovl 119921 : iKluzniak fc Rudermanl 
[1991 . 

Third, the presence of large tidal tails in which ma- 
terial will fall back onto the central object at a later 
time is a generic feature of the present set of calcula- 
tions. The mass involved in these structures is con- 
siderably larger than for binary mergers for the rea- 
sons already mentioned in Section 12.51 As the impact 
parameter in the black hole neutron star collisions in- 
creases, the mass in the tail can become even larger 
than that in the disk. Thus the properties of the fall 
back material will dominate the behavior at late times 
(|Lee. Ramirez-Ruiz fc Lopez-Ca mara 20Ql). 

The differential energy distribution in the tails is 
shown in Figure 1171 Material with negative energy is 
either in the torus surrounding the black hole, or in a 
portion of the tail that is bound to it. Fluid with positive 
energy will eventually leave the system. The thick black 
vertical line in each plot separates the torus and the 
tails morphologically, making it clear that a substantial 
portion of the tails is bound and will fall back. In 
the simple analytical estimates performed initially for 
stellar interactions with supermassive black holes, the 
distribution of mass with specific energy of the material 
from the disrupted star was constant, thus giving rise 
to a fall back accretion rate Mfh oc t~ 5 / 3 , computed 
assuming ballistic trajectories ([Reesl [19881 ) . We have 
computed the corresponding accretion rate here as well, 
also by assuming that the fluid in the tails is on ballistic 
trajectories in the potential well of the central mass (al- 
lowing us to follow it for a much longer time) . The more 
complex interaction we have outlined for the case of 
comparable masses produces a different decay law, closer 
to Mfb oc £~ 4 / 3 , see Figure [THJ There is also variability 
on shorter time scales superposed on this decay, due to 
inhomogeneities in the tidal tails. Characterizing this 
requires a full hydrodynamical and thermodynamical 
treatment of th e motion in the tails at late times 
(iRosswod 120071: ILee. Ramirez-Ruiz fc Lopez-Camaral 
2009: [Metzger et al.ll2009t r 

Essentially, the bulk of the material in the tail will 
return to the vicinity of the compact object within a few 
seconds. As it has finite angular momentum, it will not 
directly impact the primary but describe an accentric 
orbit around it. Moreover, since it is not composed of 
test particles but is a fluid stream, dissipation close to 
periastron will circularize the orbit at a radius roughly 
equal to the distance from which it was ejected in the 
first place. If angular momentum transport processes are 
present and the gas is able to cool with some efficiency, 
it will then form an accretion disk and feed the primary 
with mass and energy. The viscous time scale will be 
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Fig. 17. — Differential distribution of mass with specific energy 
for the fluid at the end of runs Lrj, Lrjl^^, Li, L2 and L2r 5 / 3 . 
The thick (thin) lines are for T = 2 (T = 5/3). The vertical line 
marks the morphological division between matter that lies in the 
torus around the black hole and that contained in the tidal tails. 
A fraction of the latter has positive energy and will escape the 
system. 
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Fig. 18. — The fall back accretion rate onto the central object 
from tidal tails (in solar masses per second) is shown for runs Lrj, 
Lor5/3, Li, L2 and L2r 5 /3, along with a reference power law with 
decay index —4/3. The thick (thin) lines are for T = 2 (r = 5/3). 

much longer than the dynamical re-injection time scale, 
and thus for transport purposes it will be as if a disk 
had been formed practically impulsively around the black 
hole (or neutron star), subsequently evolving on a secular 
time scale. It is thus possible, in principle, to account 
gcnerically in this scenario for a prompt episode of energy 
release as well as for activity many dynamical time scales 
later. 

Finally, we have shown also that the nature of the sec- 
ondary (neutron star vs. white dwarf) will lead to a 
substantially different final configuration. In the white 
dwarf case the accretion disk is much larger, and thus 
the densities are substantially lower (the total mass be- 
ing comparable). The key question in this case is whether 
the gas can find an efficient mechanism to cool, and thus 
accrete. Otherwise the material will only be blown to 
large radii and not release enough gravitational binding 
energy to account for the energy budget of a GRB. 

5.2. Time delay 

The observed event rate of SGRBs as function of 
rcdshift can provide constraints on a time delay be- 
tween t he formation of the progenitors and the ex- 
plosion dGuetta fc PiraiJ 12001 [20061 : iNakar et all [2001 
iBloom & Prochaskal 120061 1 SGRBs exhibit great di- 
versity in terms of their host galaxies, and a cur- 
sory comparison of the redshift distribution of SGRBs 
with the universal star formation rate (SFR) reveals 
what appears to be a significant time delay of a few 
Gyr. A large progenitor lifetime would help explain 
the a pparent high in c idenc e of galaxy cluster member- 
ship ([Pedersen et all 120051: IBloom et all 120061: iBergerl 
I2007T ) . while on the other hand, shorter lifetimes are 
required to explain t he population of SGRBs at mod- 
erately high redshift (jBerger et al.l 120071 : iGraham et al.l 
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Fig. 19. — The predicted cumulative redshift distribution of 
SGRBs under the assumption that they are produced by tidal cap- 
ture of two neutron stars in the cores of globular clus ters. The 
meaning of different line styles is the same as in Fig l 131 



2009). Short delays have been pointed out as a possibil- 
ity for merging bin aries directly from populations syn- 
thesis calculations (IBelczvnski. Bulik fe K alogcra 120021 : 
iBelczvnski et al.ll2006t lO'Shaughnessv et al.ll2008l ). and 
could be responsible for a significant fraction of the ob- 
served event rate. 

An alternative approach for constraining the distri- 
bution of time delays may be to use the event rates 
of SGRBs in different types of gala xies (|Gal-Yam et al.l 
I2008t iZheng fc Ramirez- Ruizl 120071) . On average, early- 
type galaxies have their stars formed earlier than late- 
type galaxies, and this difference, together with the time 
delay between progenitor formation and SGRB outburst, 
inevitably leads to different burst rates in the two types 
of galaxies. For instance, the morphological types for 
SGRBs reflect a higher incidence of early-type galaxies 
than Type la supernovae and this would suggest asso- 
ciated progenitor_l_^^ exceeding a few 
Gyrs ( Zhe ng fe Ramirez- Ruizll2007l ). Making more quan- 
titative statements about the ages of the progenitor sys- 
tems is not only hampered by small number statistics 
but also from the lack of robust predictions of the distri- 
bution of merger sites as a function of time. 

Observational evidence for a significant time delay be- 
tween tracers of star formation and SGRB outburst ex- 
cludes that the production of SGRBs is exclusively re- 
lated to the short time-scale for evolution and death of 
massive stars, which are believed to be the progenitors 
for long GRBs. However, significant time lags between 
the cosmic star formation rate and the SGRB redshift 
distribution occur naturally if the pr ogenitor is a com- 
pact binary in the fie l d that merges dCenko et al.l 120081 : 
iSalvaterra et al.ll2008t iHopman et al.ll2006D . as well as m 
the case of a close encounter between two compact stellar 
objects in a dense stellar environment. For a merger, the 
time delay reflects the time for the two compact objects 



to merge by emission of gravitational waves, while for 
close encounters a significant rate is only achieved when 
the stellar density rises significantly through the core col- 
lapse of a stellar cluster. Figure H9l shows our predicted 
cumulative redshift distribution for the progenitors of 
SGRBs assuming they are produced by tidal capture of 
two neutron stars in GCs. Further redshift determina- 
tions are required to help differentiate between various 
ways of forming a short GRB, although, as we argue next 
in Section 15.31 detailed observations of the astrophysics 
of individual GRB host galaxies may be essential before 
stringent constraints on the lifetime of short GRB pro- 
genitors can be placed. 

5.3. Galactic environment 

The observed offsets from what has been argued 
are the plausible hosts, if true, also holds important 
ramifications for the sort of viable progenitors (see 
Belczvnski et al . 2006; Zcm p. Ramirez-Ruiz fc Diemandl 



20091 ). Very large offsets seen from early-type hosts 
would seem to be at odds with progenitor systems with 
small systematic kicks such as those occurring in glob- 
ular clusters, although with such large physical offsets 
the possibility remains that the association with the pu- 
tative host is coincidental. On the other hand, based 
on the small offsets from some low-mas s gala x y hosts 
dProchaska et al.l 120061 iSoderberg et al.l 120061: iBergerl 
120091 : iBloom fc Prochaskall2006l ). SGRB progenitors can- 
not all have large systematic kicks at birth and inherently 
large delay times from formation. 

Compact binaries in the field are expected to experi- 
ence a kick, leading to mergers away from their point of 
origin. The displacement depends on the distribution of 
kick velocities, merger times, and host masses, with pre- 
dictcd values in the ra nge 10-100 kpc dFrver et al.ll"l 999: 
Belczvnski et al.ll2006t iZemp. Ramirez-Ruiz fc Diemandl 



2009f ). depending on the formation channel, host galaxy 



type and mass. By contrast, GCs are expected to show 
on average moderate displacements: the spatial distribu- 
tion of GCs peaks around the half-light radius — typically 
a few kpc — of the host galaxy. What is more, due to 
the large red giant star density in the GC core, there is 
the possibility for the interaction of the external shock 
with a denser external medium tha n that of the IGM 
(|Parsons. Ramirez-Ruiz fc Ledl2009D . 

If a significant fraction of SGRBs is indeed produced 
in GCs we also expect to see a strong preference for this 
scenario occurring in luminous host galaxies, because the 
number of GCs increas es steeply with th e host galaxy 
luminosity as oc L> 15 (jMcLaughlinl 119991 ) . whereas the 
number of primordial binaries increases only as oc Ly . 
Red, early-type galaxies are on average more luminous 
than blue, late-type galaxies, which would lead to a 
higher incidence of this kind of SGRB in ellipticals when 
compared to spiral galaxies. The current limited sample 
indicates that, if unidentified hosts are in fact ellipticals, 
SGRBs would be approximatel y evenly dist ributed be- 
tween early and late-type hosts (Bcrgcr 2009). However, 
the morphology of galaxies changes with time due to in- 
ternal (secular) evolution, and in particular when they 
merge. As a result, the fraction of late-type to early- 
type galaxies increases toward higher redshift, and hence 
SGRBs in GCs might occur also more frequently in spiral 
galaxies. 
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In other words, if (old) GCs are formed in the high- 
est d ensity peaks before reionization (e.g. iMoore et al.l 
2006), and subsequently were accreted in galaxies 
through mergers, SGRBs or any other transient con- 
nected with GCs would provide a unique tracer of the 
hierchical build-up of galaxies. This look-back on galaxy 
evolution is of course limited by the observability of such 
transients, both due to their intrinsic brightness and 
the duration of the event, as well as due to the time 
delay discussed in Section 15.21 The latter time delay 
is probably also what limits the possible contribution 
from intermediate- age GCs which are thought to form 
in gas-rich galaxy mergers. The formation time of these 
intermediate-age GCs plus the time for their cores to col- 
lapse will significantly shift the potential production of 
SGRBs towards lower redshift with respect to old GCs. 

Up to this point, we have argued for the production of 
SGRBs in GCs as oppposed to originating from primor- 
dial compact binaries within the host galaxy field, but 
we have not yet distinguised between mergers or close 
encounters between compact stellar objects in GCs. Like 
intermediate-age versus old GCs, we also expect the red- 
shift distribution of SGRBs in the merger scenario to be 
skewed toward lower redshift than in the encounter sce- 
nario. From the global encounter rate in Figure 1161 we 
expect a significant SGRB event rate only below z ~ 1. 
This is still consistent wi th SGRBs found recently around 
z ~ 1 (Be rger et al.ll2007j ) , but is very challenging or per- 
haps impossible in case of the additional time delay in 
case of merging compact binaries. 

5.4. Nuclear clusters 

The high stellar densities in which both close encoun- 
ters and three-body interactions become significant may, 
besides in the collapsed cores of GCs, also be reached 
in the nuclei of galaxies. Of particular interest are the 
nuclear (stellar) clusters (NCs) as they might also expe- 
rience core collapse. NCs have been found in 50%-75% 
of both late-type and early-type galaxies, but they are 
absent in the elliptical galaxies brighter than absolute 
magnitude Mb ~ —19 (although this could in part be 
an observational bias due to the presence of steep cen- 
tral cusps), and the frequency of nucleation also falls to 
zero in d warf galaxies fainter than Mb ~ — 12 (see the 
review bv lBokerfl2008ft NCs have typical half-light radii 
(~ 3.5 pc; B oker et al.l I2002T ) similar to GCs, but their 
typical masses (~3x 10 6 M^: IWaTcher et all 12005V and 
hence their average density is also higher than that of 
GCs by the same factor. 

On the other hand, the half-mass relaxation time of 
NCs is about an order of magnitude longer than those 
in GCs, both due to their higher typical mass as well as 
due to the higher velocity dispersion of the surrounding 
galaxy. The latter results in a flow of energy from the 
galaxy to the NC, which op poses core collapse and might 
even lead to co re expansion (lDokuchaev fc OzernoilU985l : 
lKandrutl l990) . The result is a mi nimum compactness of 
NCs in order to resist expansion (|Quinlanlll996l : iMerrittl 
2009), which is close to their observed sizes. This, to- 
gether with their longer half-mass relaxation times when 
compared with GCs, argues against NCs experiencing 
core collapse and thus providing a fertile ground for the 
production of SGRBs via close encounters between com- 
pact stellar objects. In addition, the presence of a pos- 



sible (super or in termediate) massive black halo in NCs 
(jSeth et al.l 120081 ) inhibts core collapse, and NCs with 
black holes always expand (|Merrittil2009f) . 

Still, it is possible that NCs were formed with a high 
initial concentration (and without an as yet significantly 
massive black hole) at the highest density peaks so that 
their cores collapsed early on. Later on the core and NC 
as a whole expanded, partly at least due to growing heat 
input from the surrounding galaxy of which the veloc- 
ity dispersion increases as a result of its (secular and/or 
hierarchical) evolution. An additional and likely even 
stronger expansion of a NC is expected in case the host 
(dwarf) galaxy is stripped away when it is accreted by 
a larger galaxy; a formation history that is believed to 
be a pplicable also to the m ost massive GCs in the Milky 
(e.g iGeorgiev et all 120091 ). In both cases, the currently 
observed larger size and longer half-mass relaxation time 
would not be an indication of a possible core collapse in 
the past. 

If NCs indeed experienced such a core collapse signifi- 
cantly earlier than GCs, we expect from comparing the 
earlier collapse and fiducial model in Figure [16] that the 
SGRBs potentially produced through close encounters of 
compact objects in NCs will be skewed to higher redshift 
z > 0.7. At the same time, we expect the SGRB to be lo- 
cated in the center of typically a faint (nucleated) dwarf 
galaxy. Also, the multiple st ellar populations obs erved 
in local nuclear clusters (e.g. Wal cher et al.1 12006) very 
likely imply on-going star formation at higher redshift. 
This is all in line with the recently observe d high-redshif t 
SGRBs and their host galaxy properties (|Bergerll2009j) . 
but clearly requires further study of NCs and in partic- 
ular evidence for past core collapse. 

We finally note that even though the presence of a 
central black hole inhibits core collapse and hence pre- 
vents short (two-body) relaxation times, secular resonant 
relaxation can operate on a much smaller time scale. 
This alternative relaxation process can increase the ellip- 
ticity of stellar objects on (near) Keplerian orbits, and 
hence not only bring the m closer to the black hole (e.g. 
iRauch fc Tremainel 119961 ). but also to each other. This 
increases again the change of close encounters, and, if it 
involves two compact objects, the potential to produce 
SGRBs in the vicinity of massive black holes. 

6. SUMMARY AND CONCLUSIONS 

It is evident from the work presented here that tidal 
capture in close encounters and collisions between com- 
pact stellar objects in the cores of GCs can provide a vi- 
able channel for the production of SGRBs at rates which 
are significant when compared to those from the merg- 
ers of primordial compact binaries, in both GCs and the 
field. For GCs with small primordial binary populations, 
we have argued that the formation of coalescing compact 
binaries which are capable of powering a SGRB is domi- 
nated by dynamical exchanges of NSs (BHs) with tidally 
formed WD-NS(BH) binaries and by the close encounters 
of NSs with NS/BHs. 

Through a realistic though conservative calculation, we 
predict an event rate for the tidal capture and collision 
of two neutron stars in the collapsed core of a GC that 
steeply increases to ~ 50yr _1 Gpc -3 around z ~ 0.7, 
and is followed by a gradual decline to ~ 30 yr _1 Gpc -3 
at z — 0. This is consistent with the currently observed 
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event rate and redshift distribution of SGRBs. Further- 
more, since the number of GCs both steeply increases 
with galaxy luminosity and peaks at the half-light radius 
of the host galaxy, we expect SGRBs to appear in both 
late and early-type galaxies and to be displaced from the 
galaxy center, in line with the cursory identification of 
and location in host galaxies. 

Using detailed hydrodynamics simulations we have ex- 
plored the collision of NSs with WD/NS/BHs for a 
range of system parameters, to complement previous 
knowledge on binary mergers. We have shown that 
close encounters involving neutron stars and/or black 
holes can account for both a prompt episode of en- 
ergy release, as well as for activity many dynamical 
times later through a accretion disk formed from the 
re-infall of disrupted material expelled in a bound tail. 
Relatively small changes in the impact parameter of 
the close encounter, as well as the nature of the sec- 
ondary clearly result in a significant diversity in the fi- 
nal configuration. The significant amount of material 
that is dynamically ejected comes from the multiple pe- 
riastron passages experienced before full tidal disrup- 
tion takes place, each one carrying a distinctive sig- 
nature of its thermodynamical history. Thus, while it 



could contribute signifi cantly to the observed abundance 
of r-process elements (|Lattimer fc Schrammlll974l 119761 : 
iSvmbalistv fc Schramidll982t iFreiburghaus et al.lll999ft . 
not all of it need do so in principle. 
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